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This  final  technical  report  is  composed  of  two  reprints 

f 

(attached)  plus  the— following  summary  of  those  two  papers  and 
a  discussion  of  other  work  that  is  now  in  preparation  for 
publication. 

1.  Oxygen :  It  was  noted  in  the  original  proposal  (p.  6)  that 

"we  make  corrections  to  the  fugacity  for  gas-gas  and  gas-water 
interactions  in  the  vapor  phase."  Corrections  for  molecular 
interactions  in  a  pure  gas  can  be  made  from  a  knowledge  of 
the  virial  coefficients  for  the  gas.  The  presence  of  water 
vapor,  however,  complicates  the  problem,  especially  at  higher 
temperatures,  where  the  vapor  pressure  of  water  becomes  larger. 
In  the  absence  of  experimental  knowledge  about  gas-water 
interactions,  we  were  using  theoretical  calculations  to 
estimate  the  effect  of  the  interaction.  After  it  was  found 
that  the  theoretical  methods  did  not  provide  corrections  to  the 
solubility  data  with  sufficient  accuracy,  we  finally  decided 
to  adopt  the  direct  empirical  approach  described  in  the 
attached  reprint  from  Vol.  8,  No.  9  (Benson,  Krause  and  Peterson 
1979)  of  the  Journal  of  Solution  Chemistry.  Although  additional 
time  is  required  to  carry  out  the  measurements  from  which 
the  effect  of  the  total  pressure  is  determined,  the  results 
for  the  Henry  coefficient  are  unequivocal  and  fully  corrected. 

The  precision  of  the  measurements  of  the  solubility  of 
oxygen  (random  error  less  than  0.02%)  is  an  order  of  magnitude 
better  than  previous  measurements  on  any  gas  over  a  significant 
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temperature  range,  and  tests  and  comparison  with  other  work 
indicate  that  systematic  errors  probably  are  negligible. 
Furthermore  the  accuracy  of  the  measurements  has  made  it 
possible  to  show  that  the  expression 

In  k  =  aQ  +  a^/T  +  a2/T2 

is  a  much  better  way  to  express  the  variation  of  solubility 

with  temperature  than  the  relationships  previously  used. 

Mass  spectrometric  measurements  of  the  isotopic  fractionation 
34  32 

of  02/  02  during  solution  have  provided  the  first 

quantitative  determination  of  the  variation  of  the  fraction¬ 
ation  with  temperature,  and  have  made  possible  an  estimation 
of  the  size  of  the  cavity  occupied  by  the  oxygen  molecule 
in  the  water. 

A  second  manuscript,  "The  concentration  and  isotopic 
fractionation  of  gases  dissolved  in  freshwater  in  equilibrium 
with  the  atmosphere:  1.  Oxygen,"  has  been  published  in 
"Limnology  and  Oceanography"  (Benson  and  Krause,  1980).  A 
copy  is  included  here.  In  it  we  have  discussed  how  to 
obtain  accurate  values  for  the  concentration  of  an  atmospheric 
gas  in  water  in  equilibrium  with  air,  including  corrections 
for  molecular  interactions  in  the  vapor  phase.  The  results 
for  oxygen  show  that  the  UNESCO  tables  for  the  solubility 
of  oxygen  give  values  in  freshwater  that  are  low  by  from  0.157, 
to  0.35%  depending  upon  temperature. 
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2 .  Hydrogen :  We  have  carried  out  several  sets  of  determina¬ 

tions  of  the  solubility  of  hydrogen.  The  results  give  very 
accurate  values  for  oceanographic  and  limnological  work,  and  they 
are  especially  significant  because  of  the  information  they  provide 
about  the  properties  of  liquid  water. 

A)  Pure  Water:  Fig.  1  shows  the  results  of  our 

measurements  on  hydrogen  in  distilled  water  from  0  to 
60°C.  The  clusters  of  solid  circles  at  0,  15,  35  and 
55°C  are  drawn  to  indicate  that  measurements  at  several 
pressures  were  made ,  from  which  the  corrections  for 
molecular  interactions  in  the  vapor  phase  were  obtained. 

(All  points  really  lie  essentially  on  the  line  at  the 
same  temperature.)  Although  the  virial  coefficients  for 
pure  hydrogen  show  that  it  is  a  "repulsive  gas,"  the 
hydrogen-water  vapor  system  behaves  very  nearly  like  an 
ideal  gas.  In  fact,  the  variation  of  \  with  t  (°C)  is 
given  by 

\  =  0.000345  -  1.83  x  10_5t, 

which  shows  that  below  approximately  room  temperature  it 
is  a  slightly  "attractive  gas,"  but  it  becomes  "repulsive" 
above  room  temperature.  Comparison  with  Eq .  21  in  Benson, 
Krause  and  Peterson  (1979)  shows  that  relative  to  oxygen 
the  vapor  phase  correction  is  smaller  and  varies  less 
rapidly  with  temperature.  Again,  our  function 
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provides  an  excellent  fit  for  the  data.  The  standard 
deviation  for  all  the  pure  water  measurements  from  0 
to  60°C  is  0.016%.  Fig.  1  also  shows  the  experimental 
results  of  Crozier  and  Yamamoto  (1974)  and  Gordon  et  al. 
(1977).  The  large  scatter  of  their  data  is  evident, 
together  with  systematic  deviations  from  our  results. 

Note  that  the  ordinate  axis  has  been  compressed  by  a 
factor  of  two  compared  to  Fig.  4  in  BKP  for  oxygen,  in 
order  to  keep  the  hydrogen  points  on  the  graph.  It  is 
possible  to  calculate  the  thermodynamic  function  changes 
for  the  solution  of  hydrogen  in  water,  as  we  have  done 
for  oxygen,  but  these  will  not  be  discussed  here. 

B)  Seawater :  Although  it  is  well  known  that  gas  solubilities 

decrease  with  increasing  salinity  (the  "salting-out"  effect)  , 
the  functional  dependence  of  solubility  on  S  has  been 
difficult  to  determine  because  very  accurate  measurements 
arc  required.  Carpenter  (1966)  found  a  non-linear 
variation  for  the  solubility  of  02  with  chlorinity,  and 
he  used  a  quadratic  relationship.  Green  and  Carritt  (1967) 
reported  a  series  of  measurements  on  oxygen  to  obtain  the 
r.rio  of  the  solubility  in  seawater  to  that  in  pure  water. 
The  design  of  the  experiment  made  potential  systematic 
err<  in  the  individual  solubility  determinations  tend 
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within  the  limits  of  their  measurements  the  Setschenow 
relationship  correctly  describes  the  dependence  of  the 
oxygen  ratio  on  chlorinity.  Furthermore,  their  results 
provided  the  first  reasonably  accurate  determination  of 
the  way  the  salting-out  coefficient  varies  with  temperature. 

We  have  measured  the  dependence  of  solubility  on 
salinity  up  to  approximately  50%*  at  0 ,  10,  20  and  35°C. 

The  results  are  shown  in  Fig.  2  where  it  is  clear  that 
In  (kg/ko)  is  extremely  linear  with  salinity.  (Fits  to 
salinity  and  chlorinity  are  equally  good.)  Although  the 
straight  lines  in  Fig.  2  suggest  that  the  Setschenow 
equation  fits  the  results  very  well,  the  tightness  of  fit 
does  not  really  become  apparent  until  it  is  recognized 
that  the  precision  of  the  measurements  is  such  that  even 
the  slopes  of  the  four  lines  are  known  with  great  accuracy 
(0°C :  0.46%,  10°C:  0.12%,  20°C:  0.15%,  35°C:  0.04%). 

In  other  words ,  these  very  accurate  new  results  for 
hydrogen  not  only  verify  the  applicability  of  the 
Setschenow  relationship  to  hydrogen  in  seawater,  but  also 
now  make  it  possible  to  determine  the  functional  dependence 
of  the  salting-out  coefficient  upon  temperature.  This 
is  illustrated  in  Fig.  3  where  the  salting-out  coefficient 
is  plotted  vs.  1/T.  (Error  bars  are  shown  on  the  solid 
circle  at  0°C,  but  the  errors  for  the  other  three  poix*ts 
are  less  than  the  size  of  the  circles.)  It  is  clear  that 
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the  points  lie  on  a  smooth  curve,  and  -  like  In  k  -  it  is 
parabolic  in  1/T.  The  open  squares  representing  the 
results  of  Gordon  et  al .  are  generally  consistent  with  our 
values,  although,  of  course,  their  much  larger  errors  and 
scatter  preclude  determination  of  the  temperature  dependence 
with  any  certainty. 

Our  final  equation  for  the  variation  of  the  Henry 
coefficient  with  both  temperature  and  salinity  is 

In  k  =  4.13352  +  4672. 47/T  -  768235/T2 

+  S(0. 026112  -  15.684/T  +  2737. 0/T2). 

The  salting-out  coefficient  in  the  parenthesis  has  its 
minimum  at  a  temperature  not  very  different  from  that  for 
the  solubility  in  pure  water.  The  equation  used  by  Green 
and  Carritt  for  the  Bunsen  coefficient  of  oxygen  was  of  the 
form 


In  3  =  A  +  B/T  +  C  In  T  +  DT 

+  Ci< A'  +  B'/T  +  C'lnT  +  D’T). 

Gordon  et  al .  fitted  their  hydrogen  data  together  with  those 
of  Crozier  and  Yamamoto,  to  the  following  function 

In  8  =  A  +  B/T  +  C  In  T 

+  S (A '  4-  B'T  +  C'T2) . 

The  fit  of  all  our  data  (for  both  pure  water  and  seawater) 
to  our  function  is  shown  in  Fig.  4.  There  were,  of  course, 
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many  more  points  at  zero  salinity  than  those  indicated. 

The  values  of  Gordon  et  al .  also  are  shown  as  a  function 
of  salinity  at  five  temperatures  from  0  to  29°C.  In 
addition  to  large  fluctuations  they  exhibit  apparent 
systematic  deviations  which  depend  upon  temperature. 

Fig.  5  compares  the  experimental  values  of  both  Gordon  et 
al.  and  Crozier  and  Yamamoto  with  our  results  for  the 
temperature  variation  of  the  solubility  of  hydrogen  in  water 
with  the  30-407,,  salinity  range  typical  of  the  oceans. 

3.  Data  accumulation  and  analysis  system. 

For  very  precise  isotopic  ratio  measurements  with  the 
mass  spectrometer  we  have  designed  and  constructed  a  very  con¬ 
venient  new  system  for  taking  and  analyzing  the  data.  It 
incorporates  a  voltage-to- frequency  converter  and  commercially 
available  timer,  counter,  and  printer  units  into  an  automatic 
sample  changing  system  that  electronically  averages  noise 
from  the  difference  amplifier  in  the  null  balance  system. 

The  system  yields  superb  results.  The  precision  achieved 
with  electronic  averaging  of  the  noise  is  both  greater  and 
more  objective  than  with  visual  averaging.  In  addition,  the 
time  for  taking  and  analyzing  the  data  is  reduced  by  approxi¬ 
mately  307,. 
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A  very  precise  and  accurate  new  method  is  described  for  determination  of  the 
Henry  coefficient  k  and  the  isotopic  fractionation  of  gases  dissolved  in  liquids. 
It  yields  fully  corrected  values  for  k  at  essentially  infinite  dilution.  For  oxygen 
the  random  error  for  k  is  less  than  0.02%,  which  is  an  order  of  magnitude 
better  than  the  best  previous  measurements  on  that  or  any  other  gas.  Extensive 
tests  and  comparison  with  other  work  indicate  that  systematic  errors  probably 
are  negligible  and  that  the  accuracy  is  determined  by  the  precision  of  the 
measurements.  In  the  virial  correction  factor  ( I  +  A P, ) .  where  P,  is  the  total 
pressure  of  the  vapor  phase,  the  coefficient  A  for  oxygen  empirically  is  a  linear 
decreasing  function  of  the  temperature  over  the  range  0-60  C.  The  simple 
three-term  power  series  in  I  IT  proposed  by  Benson  and  Krause, 

In  k  =  ao  +  aJT  +  afiT2 

provides  a  much  better  form  for  the  variation  of  k  with  temperature  than  any 
previous  expression.  With  a0  =  3.71814,  a,  -  5596.17,  and  a2  —  -1049668, 
the  precision  of  fit  to  it  of  37  data  points  for  oxygen  from  0-60  C  is  0.018% 
(one  standard  deviation).  The  three-term  series  in  1IT  also  yields  the  best  fit 
for  the  most  accurate  data  on  equilibrium  constants  for  other  types  of  systems, 
which  suggests  that  the  function  may  have  broader  applications.  The  oxygen 
results  support  the  idea  that  when  the  function  is  rewritten  as 

Ink  =  -(A,  +  A2)  +  +  A2(^Y 

it  becomes  a  universal  solubility  equation  in  the  sense  that  A  2  is  common  to  all 
gases,  w  ith  Tx  ami  .4,  characteristic  of  the  specific  gas.  Accurate  values  are 
presented  for  the  partial  molal  thermodynamic  function  changes  for  the 
solution  of  oxygen  in  water  between  the  usual  standard  states  for  the  liquid 
and  vapor  phases.  These  include  the  change  in  heat  capacity ,  which  varies 
inversely  with  the  square  of  the  absolute  temperature  and  for  which  the  random 
error  is  0.15 To-  Analysis  of  the  high-temperature  data  of  Stephan  et  a!.%  in 
combination  with  our  values  from  0-60° C,  shows  that  for  oxygen  the  four- 
term  series  in  1  IT, 

)n  k  =  -4.1741  +  J.3104  x  10 NT  -  3.4170  x  \0«ITJ  +  2.4749  x  10 VT* 
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where  p  —  kx  and  p  is  the  partial  pressure  in  atmospheres  of  the  gas,  probably 
provides  the  best  and  easiest  way  presently  available  to  calculate  values  for  k 
in  the  range  100-288  C,  but  more  precise  measurements  at  elevated  tem¬ 
peratures  are  needed.  The  new  method  permits  direct  mass  spectrometric 
comparison  of  the  isotopic  ratio  O p  ,  in  the  dissolved  gas  to  that  in  the 
gas  above  the  solution.  The  fractionation  factor  a  =  33k/3*k  varies  from 
approximately  / .00085  (  ±0.00002)  at  0  C  to  / .00055  (  ±  0.00002 )  at  60  C. 
Although  the  results  provide  the  first  quantitative  determination  of  a  vs. 
temperature  for  oxygen,  it  is  not  possible  from  these  data  to  choose  among 
several  functions  for  the  variation  of  In  a  with  temperature.  If  the  isotopic 
fractionation  is  assumed  to  be  due  to  a  difference  in  the  zero-point  energy  of 
the  two  species  of  oxygen  molecules,  the  size  of  the  solvent  cage  is  calculated 
to  be  approximately  2.5  A.  The  isotopic  measurements  indicate  that  sub¬ 
stitution  of  a  "O ,  molecule  for  a  J102  molecule  in  solution  involves  a  change 
in  enthalpy  with  a  relatively  small  change  in  entropy. 


KEY  WORDS:  Gas  solubility:  isotopic  fractionation,  oxygen;  water; 
temperature  dependence  ;  Henry  coefficient;  thermodynamic  functions. 


1.  INTRODUCTION 

Although  the  solution  of  gases  in  water  has  been  studied  for  approximately 
two  centuries,  obtaining  accurate  values  for  gas  solubilities  has  proved  to  be 
difficult.  The  history  of  the  subject  has  been  reviewed  thoroughly  by  Markham 
and  Kobe,11’  Battino  and  Clever.12'  Wilhelm,  Battino,  and  Wilcock,131  and 
Franks.14'  In  this  paper  the  discussion  of  earlier  work  on  oxygen  will  be  brief, 
except  for  a  comparison  of  the  results  presented  here  with  the  best  previous 
values. 

Because  a  knowledge  of  its  solubility  is  important  for  many  kinds  of 
work,  oxygen  has  been  studied  more  extensively  and  carefully  than  any  other 
gas.  In  retrospect,  the  early  measurements  by  Winkler15,6'  and  Fox17’  during 
the  last  two  decades  of  the  19th  century  and  the  first  decade  of  the  20th  must 
be  considered  very  good  for  that  period.  They  survived  as  standard  values  in 
the  Handbook  of  Chemistry  and  Physics**'  (Winkler),  as  standards  for  oceano¬ 
graphic  work  (Fox),  and  as  standards  for  sanitation  engineering  (the  oxygen 
values  of  Fox  recalculated  by  Whipple  and  Whipple19')  for  a  half  century 
and  longer.  The  accuracy  of  the  Fox  and  Winkler  values  was  brought  into 
question  in  1955  by  the  oxygen  determinations  of  Truesdale,  Downing,  and 
Lowden.110’  Morrison  and  Billett,111’  Steen,112’  Morris,  Stumm,  and  Galal,113’ 
and  Elmore  and  Hayes1141  carried  out  new  analyses  on  oxygen,  and  in  1961 
Benson  and  Parker1151  used  mass  spectrometric  measurements  of  the  solu¬ 
bility  ratios  Na/Ar  and  Na/02  to  interrelate  and  evaluate  existing  data.  The 
differences  among  workers  indicated  that  the  uncertainties  in  the  solubility 
values  were  approximately  1-2%. 

The  need  for  better  solubility  values,  primarily  for  the  study  of  oxygen 
dissolved  in  natural  waters,  generated  a  series  of  attempts116-22’  to  achieve 
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higher  accuracy.  By  1970,  as  we  shall  see  later,  four  laboratories,  using 
different  methods  of  equilibration  and  both  gasometric  and  Winkler  analytical 
techniques,  had  obtained  agreement  over  the  range  0-35  C  within  the  scatter 
of  their  data,  i.e..  approximately  ±0.2°^.  Together,  these  constitute  the  most 
accurate  solubility  data  currently  available  for  any  gas. 

Various  methods  have  been  used  to  express  solubility  values  as  a  function 
of  temperature.  Most  recently,  Wilhelm,  Battino,  and  Wilcock,131  following 
Clarke  and  Glew,'23’  have  fitted  a  selected  group  of  data  for  oxygen  (and 
other  gases)  to  a  four-term  series  of  the  form 

In  k  =  a0  +  flj/T  +  a2  In  T  +  a3T  H -  (I) 

where  k  is  the  Henry  coefficient.  This  series  is  essentially  an  extension  of  the 
three-term  function  first  derived  by  Valentiner.'24’  Weiss'25’  fitted  the  Bunsen 
coefficient  data  of  Carpenter'20’  and  Murray  and  Riley,122’  for  both  pure 
water  and  seawater,  to  a  three-term  series  in  the  temperature  of  the  Valentiner 
form,  with  additional  terms  for  the  effect  of  salinity.  The  Weiss  values  are 
often  used  for  oceanographic  and  other  purposes. 

During  the  development  of  the  new  apparatus  and  procedures  to  be 
described  below,  many  preliminary  measurements  were  carried  out  on 
several  gases  to  test  for  random  and,  especially,  systematic  errors.  Although 
we  now  know  there  was  a  systematic  error  in  those  preliminary  determinations 
(see  Appenfiix  A),  and  although  they  were  crude  compared  to  the  measure¬ 
ments  now  possible,  analyses  of  the  data  from  them,  and  from  Murray. 
Riley,  and  Wilson.126-271  led  Benson  and  Krause'28’  to  propose  series  of  the 
form 

In  k  =  a0  +  a, IT  +  a2[T2  +  a3/T3  +  ■■■  (2) 

for  the  variation  of  k  with  temperature.  Himmelblau'29’  had  considered  a 
function  of  this  form,  but  discarded  it  in  favor  of  a  general  six-term  equation 
of  the  second  degree  in  In  k  and  l/T.  Benson  and  Krause  found  that  a  simple 
three-term  series  in  I /T  not  only  fitted  the  preliminary  data  better  than  other 
expressions,  but  also  led  to  the  recognition  of  probable  relationships  between 
the  molecular  properties  of  the  pure  gases  and  the  thermodynamic  properties 
of  their  solutions. 

Relatively  little  is  known  about  the  isotopic  fractionation  of  gases  during 
solution.  It  was  first  observed  by  Klots  and  Benson  ,3t”  in  studies  of  oxygen 
and  nitrogen.  For  each  gas  it  was  found  that  the  heavier  molecular  species 
(3402  or  29N2)  was  more  soluble  than  the  respective  lighter  one  (3202  or 
28N2)  by  approximately  0.8  per  mil  at  OX.  and  the  effect  decreased  to  approxi¬ 
mately  0.55  per  mil  at  27X.  Polgar1311  obtained  results  for  argon  that  were 
roughly  comparable,  and  Weiss132’  observed  the  larger  effect  (approximately 
10  per  mil)  to  be  expected  for  helium.  Kroopnick  and  Craig'33’  found  an 
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effect  with  oxygen  dissolved  in  seawater  that  was  similar  to  the  results  of 
Klots  and  Benson  in  pure  water.  None  of  these  measurements  was  precise 
enough  to  provide  quantitative  information  about  the  variation  of  the 
isotopic  fractionation  with  temperature  and  salinity.  Muccitelli  and  Wen'311 
measured  the  solubility  of  D2  in  pure  water.  By  combining  their  results  with 
those  of  Crozier  and  Yamamoto'351  for  H*  in  pure  water,  they  found  the 
fractionation  ranged  from  approximately  86  per  mil  at  0°C  to  65  per  mil  at 
30X.  To  carry  out  these  calculations  they  first  fitted  each  set  of  data  separately 
to  the  generalized  Valentiner  equation  [Eq.  (I)]  to  obtain  smooth  curves. 

The  work  reported  here  is  the  first  in  a  series  designed  to  produce  very 
accurate  values  for  the  solubility  and  isotopic  fractionation  of  gases  in  pure 
water  from  0  to  60°C.  These  will  reveal  the  dependence  of  solubility  on 
temperature,  yield  information  for  a  better  understanding  of  molecular 
interactions  in  liquid  water,  and  provide  reference  values  for  geochemical 
work.  We  have  chosen  to  study  oxygen  first,  to  permit  comparison  of  results 
from  the  new  method  with  the  best  previous  determinations. 

2.  EXPERIMENTAL  METHODS  AND  THEORY 
2.1.  The  Henry  Coefficient 

The  essential  features  of  any  apparatus  for  determination  of  gas  solu¬ 
bility  are  provision  for  equilibration  of  the  desired  gas  and  liquid,  and 
means  for  finding  the  fugacity  of  the  gas  in  the  vapor  phase  and  the  mole 
fraction  of  the  dissolved  gas  in  the  liquid  phase.  We  have  used  several  different 
equilibration  arrangements.  The  best,  and  the  one  used  for  the  measurements 
reported  here,  is  shown  in  Fig.  1.  A  brief  description  of  the  equilibrator  used 
in  the  preliminary  work 1281  is  given  in  Appendix  A,  with  a  discussion  of 
evidence  suggesting  that  potential  problems  are  associated  with  methods  of 
equilibration  that  involve  thin  films  of  liquid. 

The  equilibration  system  in  Fig.  1  consists  of  four  parts,  the  spherical 
equilibrator,  the  centrifugal  pump,  the  vapor-phase  sample  bulb,  and  the 
liquid-phase  sample  bulb.  The  vapor-phase  sample  bulb  is  attached  by  means 
of  a  14/35  standard  taper  joint.  Silicone  vacuum  grease  is  used  to  seal  both 
the  joint  and  the  4-mm  pressure-type  stopcock  on  the  bulb.  The  other  joints 
utilize  standard  O-ring  compression,  with  O-rings  made  of  ethylene-propylene 
rubber.  The  equilibrator  is  composed  of  concentric  spheres  with  inner  and 
outer  volumes  of  1  and  2  liters,  respectively.  The  reentrant  connection  to  the 
inside  of  the  inner  sphere  provides  a  compact  way  of  including  a  relatively 
large  volume  for  the  vapor  phase,  so  that  dissolution  reduces  the  pressure  less, 
especially  with  the  more  soluble  gases.  For  very  expensive  gases  this  feature 
is  deleted.  The  pump  employs  a  Teflon-covered  magnetic  stirring  bar  as  an 
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Fig.  1.  Glass  apparatus  for  equilibration  or  gases  with  liquids,  with  spherical  equilibrator, 
centrifugal  pump,  and  vapor-phase  and  liquid-phase  sample  bulbs.  Arrows  indicate  the 
direction  of  liquid  flow. 


impeller.  The  stopcocks  on  the  liquid-phase  sample  bulb  are  of  the  Teflon 
plug  type  (Ace  #8194),  but  with  the  Teflon  plugs  replaced  by  special  ones  of 
either  Type  304  stainless  steel  or  glass,  with  O-ring  seals  designed  to  make 
possible  isolation  of  a  very  accurately  known  and  reproducible  volume  in  the 
bulb  when  the  stopcocks  are  closed.  For  measurements  below  room  tem¬ 
perature,  a  small  O-ring  sealed  piston-cylinder  is  set  into  the  wall  of  the  bulb 
above  the  lower  stopcock.  When  released,  after  volumetric  sampling,  this 
permits  expansion  of  the  solution  to  occur  when  the  bulb  is  removed  from  the 
bath.  The  whole  assembly  is  clamped  in  a  frame  of  welded  stainless  steel 
tubing. 

Distilled  water  from  a  central  stainless  steel  still  is  redistilled  in  a  glass 
system  (Corning  Model  AG- lb).  The  resulting  water  has  a  resistivity  greater 
than  1.5  M£2-cm.  Variations  in  the  isotopic  constitution  of  natural  waters 
are  too  small  to  have  any  detectable  effect  on  solubility  measurements,  even 
at  the  level  of  precision  reported  here,  and  no  attempt  is  made  to  analyze  the 
water  isotopically.  After  the  second  distillation,  the  water  is  degassed  in  a 
2-liter  flask  to  less  than  1  ppm  of  the  air  that  would  be  dissolved  in  it  at  room 
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temperature,  by  stirring  with  a  glass-covered  magnetic  stirring  bar  and 
periodically  pumping  on  the  flask  through  a  cold  trap.  Then  with  the  flask  of 
degassed  water  attached  to  the  filling  port  of  the  equilibrator  by  means  of  a 
"Y"  coupling,  and  after  the  coupling  and  the  equilibrator  are  carefully 
evacuated  (and  the  Teflon-covered  magnetic  stirring  bar  in  the  pump 
degassed),  the  water  is  allowed  to  flow  down  into  the  equilibrator  and  the 
liquid-phase  sample  bulb  until  the  pump  is  covered  and  the  water  level  is 
high  enough  in  the  equilibrator  so  that  dripping  will  not  occur  later  during 
pumping.  Research  grade  gas  (the  oxygen  supplied  by  Linde  was  stated  to 
have  a  total  impurity  less  than  40  ppm)  is  introduced  through  a  dry  ice- 
acetone  cooled  trap  and  the  “Y”  coupling  in  such  a  way  as  to  minimize 
backstreaming  of  water  vapor  toward  the  tank  of  gas. 

After  closing  the  filling  port,  the  assembly  is  entirely  immersed  in  a  large 
constant-temperature  bath.  The  speed  of  the  centrifugal  pump  is  adjusted 
so  that  water  flows  up  the  central  tube  and  smoothly  out  over  the  inner 
sphere  without  turbulence.  After  the  desired  time  interval,  a  known  volume 
of  the  gas  above  the  solution  is  isolated  in  the  vapor-phase  sample  bulb,  and 
a  known  volume  of  the  solution  is  closed  off  in  the  liquid-phase  sample 
bulb. 

While  extraction  of  the  gas  dissolved  in  the  water  sample  is  being  carried 
out,  the  undissolved  gas  sample  is  dried  in  dry  ice-acetone  and  liquid- 
nitrogen  cooled  traps  and  Toepler-pumped  into  a  mercury  manometric 
system  designed  for  this  work.  The  system,  including  the  Toepler  pump,  is  in 
a  thermostatic  enclosure  that  permits  access  to  the  manometer  but  maintains 
temperature  uniformity  to  0.01  °C  and  constancy  to  a  few  millidegrees  over 
several  hours.  A  special  meniscus  lighting  system  makes  it  possible  to  deter¬ 
mine  the  position  and  height  of  the  menisci  with  an  accuracy  of  0.01  mm 
using  a  Beck  Model  1000  cathetometer.  The  volume  of  the  manometer  is 
easily  changed  to  provide  the  best  compromise  between  volume  and  pressure 
for  maximum  precision.  After  measurement,  a  sample  of  the  gas  may  be 
transferred  to  a  bulb  for  analysis  with  the  mass  spectrometer. 

To  remove  the  dissolved  gas,  the  liquid-phase  sample  bulb  is  attached  to 
one  arm  of  a  2-liter  extraction  flask  containing  a  glass-covered  magnetic 
stirring  bar.  The  other  arm  is  connected  to  a  dry  ice-acetone  trap,  a  liquid- 
nitrogen  trap,  and  a  Toepler  pump,  in  succession,  with  stopcocks  between 
each  section.  A  roughing  pump  and  a  diffusion  pump,  each  trapped  for  oil 
vapor,  and  a  thermocouple  vacuum  gauge  reading  directly  to  2  x  10~4  torr 
are  also  connected  to  the  manifold.  With  the  whole  manifold  fully  evacuated, 
the  water  in  the  bulb  is  allowed  to  flow  downward  into  the  extraction  flask. 
A  cross-link  between  the  two  arms  of  the  flask  makes  it  possible  to  Toepler- 
pump  the  residual  gas  in  the  liquid-phase  sample  bulb  through  the  drying 
traps  into  a  gas  storage  bulb.  Then  “slugs"  of  gas  are  taken  from  the  vigor- 
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ously  stirred  extraction  flask,  dried  successively  in  the  two  traps,  and  then 
Toepler-pumped  into  the  storage  bulb.  After  extraction  is  complete,  the  gas  is 
transferred  to  the  manometer  and  the  number  of  moles  is  measured. 

Henry’s  law  for  a  two-component  system  composed  of  a  pure  gas  and 
water  states  that  at  constant  temperature  the  fugacity  of  the  gas  in  the  vap  -r 
phase  is  proportional  to  the  mole  fraction  of  the  gas  in  solution,  i.e., 


/  = 


k 

n  w  +  na 


(3) 


where  k  is  the  Henry  coefficient,  and  nd  and  nw  are  the  numbers  of  moles  of 
dissolved  gas  and  water,  respectively,  in  a  volume  K,  of  solution,  i.e.,  the 
volume  of  the  liquid-phase  sample  bulb.  From  the  manometric  measure¬ 
ments,  Pd ,  Vd,  and  Td,  on  the  dissolved  gas,  we  obtain 


„  P*V* 

A  RTdZd 


(4) 


where  Zd  is  the  compressibility  factor  for  the  pure  gas  at  Pd  and  Td. 

The  volume  of  the  solution.  Ks,  is  equal  to  the  volume  of  the  pure  water 
in  it  plus  the  partial  molal  volume  of  the  gas  times  nd,  i.e.. 


,,  n.M 

F5  =  -  +  l'dni 

P 


where  M  is  the  molecular  weight  of  water  and  p  is  its  density  at  the  equilibrium 
temperature  T.  With  algebraic  rearrangement, 

nd  =  M  PdVd 
n*  +  na  pV,  RTdZd 

where 

V,  =  Ts(l  -  {(-•„  -  M I p)(nd! V\)} 

Determination  of  the  fugacity  is  slightly  more  involved  For  the  moment 
we  shall  assume  that  the  Lewis  and  Randall  rule  may  be  applied,  i.e.. 


(5) 

(6) 


(7) 


where  nd  is  the  number  of  moles  of  undissolved  gas  in  the  volume  T„  of  the 
vapor-phase  sample  bulb,  and  n,  is  the  total  number  of  moles  of  gas  and  water 
vapor  in  that  volume.  The  quantity/*  is  the  fugacity  that  the  pure  gas  would 
have  at  the  equilibrium  temperature  if  its  pressure  were  the  same  as  the  total 
pressure  PK  of  the  gas  and  water  vapor.  Therefore. 


Px  exp[  -  epx\ 
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because  for  low  pressures  Z  =  I  -  9P,  where  9  is  determined  by  the  second 
virial  coefficient  of  the  pure  gas.  If  the  exponential  is  expanded  as  a  power 
series  with  higher-order  terms  neglected,  we  have 

f*  =  PA  I  -  »Pi)  =  P,Z*  (8) 

where  Z*  would  be  the  compressibility  factor  for  the  pure  gas  at  the 
equilibrium  temperature  and  the  total  pressure  Pt.  If  Z,  is  the  compressibility 
factor  for  the  vapor  phase. 


Pi  = 


nlRTZl 

V* 


(9) 


Furthermore,  from  the  nianometric  measurements  on  the  undissolved  gas. 


Substituting  this  and  tqs.  (8)  and  (9)  into  Eq.  (7)  yields 


/  = 


pay  u  t 
T,,Z,, 


Z,Z’ 


(10) 


We  may  write 

ztz*  =  (i  -  0t/\)(  i  -  «\) 

where  9,  is  determined  by  the  second  virial  coefficient  of  the  gas-water  vapor 
mixture  If  the  term  9ldP?  is  neglected,  and  we  define  O'  =  0t  +  0 ,  we  have 


and 


Z,Z*  =  1  -  6'Pt 


PVT 

f~TTV-U-  9'pJ 

1  u^.u  v  b 


Substituting  Eqs.  (5)  and  (12)  into  Eq.  (3)  yields 


k  = 


R_r\RK  T 

M  nAR  V b  9 


(I  -  B'Pt) 


or 


A-  = 


J^KKTtZtK  oT 

M  P„Va  T UZU  V b  9 


(1  -  O' Pi) 


(11) 

(12) 


(13a) 


(13b) 


If  the  molecules  in  the  vapor  phase  did  not  interact  with  each  other,  i.e.,  if  the 
vapor  phase  behaved  like  an  ideal  gas.  then  Zt  =  Z*  =  1  and  9X  =  9'  =  0, 
and  the  measured  quantity  in  brackets  in  Eqs.  (13a)  and  (13b)  would  be  a 
constant  equal  to  the  Henry  coefficient.  With  molecular  attraction  occurring 
in  the  vapor  phase,  however,  one  would  expect  the  ratio  naR[naR  to  be  larger 
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than  for  the  ideal-gas  situation  and,  correspondingly,  the  factor  (1  —  d'P t) 
would  be  less  than  unity.  Nevertheless,  if  accurate  values  were  known  for  both 
6  and  6t,  Eq.  (13)  could  be  used  to  calculate  k  from  the  measured  quantities 
including  Pv 

For  the  work  reported  here,  a  more  direct  empirical  approach  has  proved 
to  be  more  fruitful.  Suppose  Eq.  (13a)  is  rewritten  as 

£SsrH-‘"  +  "'>  "4) 

where,  again,  the  fact  that  6'Pl  «  1  has  been  used.  Equation  (14)  predicts 
that  at  a  fixed  temperature  the  quantity  in  brackets  should  be  a  linear 
function  of  Pt ,  with  k  the  intercept  when  Pl  is  extrapolated  to  zero.  By  making 
such  a  plot  versus  pressure  at  several  temperatures,  the  variation  of  6'  with 
temperature  may  be  found  empirically.  Then  one  set  of  measured  quantities, 
including  Pt,  permits  calculation  of  k. 

It  is  interesting  to  note  that  Eq.  (14)  can  be  obtained  without  invoking 
the  Lewis  and  Randall  rule  [Eq.  (7)]  and  without  carrying  the  real-gas 
corrections  throughout  the  analysis,  except  in  the  manometry.  Suppose  it  is 
first  assumed  that  the  vapor  behaves  like  an  ideal  gas.  This  leads  to  an 
equation  of  the  same  form  as  Eq.  (14),  but  without  the  factor  (1  +  8'PX). 
Then  recognizing  that  the  real-gas  behavior  of  the  vapor  phase  makes  the 
quantity  in  brackets  depend  slightly  upon  Pt,  that  quantity  can  be  expanded 
as  a  Taylor  series  in  Pt.  Equation  (14)  is  the  result  if  terms  in  higher  powers  of 
Pl  are  negligible. 

This  method  for  solubility  determinations  has  many  good  features.  In 
the  first  place,  it  yields  fully  corrected  values  for  k  at  essentially  infinite 
dilution,  the  ideal  condition  for  the  validity  of  Henry's  law.  The  method  is 
straightforward  and  involves  measurements  that  can  be  made  with  both  very 
high  precision  and  accuracy.  Negligible  error  is  introduced  into  k  if  the 
approximate  expression 

Px  =  +  f.v  (15) 

is  used  for  Pt  in  the  small  correction  term.  The  quantity  PKV  is  the  saturated 
vapor  pressure  of  water  at  the  temperature  of  equilibrium,  and  Z0  is  the 
compressibility  factor  for  the  pure  gas  at  its  equilibrium  partial  pressure. 
Equation  (15)  assumes  that  the  gas  and  water  molecules  do  not  interact  in 
the  vapor  phase.  No  other  corrections  for  water  vapor  are  required,  because 
all  manometric  measurements  are  made  on  dry  gases. 

It  is  not  necessary  to  have  accurate  values  for  the  partial  molal  volume  of 
gases  with  low  solubilities  because  the  difference  between  V't  and  is  very 
small.  With  oxygen,  for  example,  the  largest  difference  is  approximately 


664 


Benson.  Krause,  and  Peterson 


30  ppm  (at  0  C  and  1  atm).  In  fact,  with  the  empirical  approach  discussed  in 
relation  to  Eq.  (14),  negligible  error  is  introduced  by  substituting  T,  for  K'. 
This  may  be  seen  from  Eq.  (6):  F,  =  K,[  1  -  (fd  —  Mjp){nAjV„)\  is  approxi¬ 
mately  equal  to  K,[l  -  (rd  -  M/p)fi(P^  -  /»„)],  where  /9  is  the  Bunsen 
coefficient  and  /? (r'd  -  M/p)  =  y  is  a  constant  at  a  given  temperature.  After 
substitution  for  V,  in  Eq.  (14)  and  rearrangement,  one  obtains 


JL  'hjl  ij?  y 

4/  ntR  yh  P 


*1(1  -  y^wv)  +  (y  +  <?  )/">) 


(16) 


Therefore,  when  the  quantity  in  brackets  on  the  left  is  plotted  versus  part 
of  the  V'a  correction  is  subsumed  in  the  empirical  determination  of  the  factor 
multiplying  Pt.  Assuming  that  yP„y  is  negligible  compared  to  unity  (for 
oxygen  from  0  to  60  C,  the  largest  error  introduced  is  only  3  ppm,  but  even 
for  a  more  soluble  gas  like  xenon  the  error  would  be  less  than  50  ppm),  Eq. 
(16)  becomes 


k' 


M  ndR  l\pl 


=  k(l  +  A/\) 


(17) 


where  A  =  y  +  8' . 

The  design  of  the  experiment  minimizes  the  significance  of  almost  all 
systematic  analytical  errors.  In  Eq.  (13b)  note  that  the  eight  manometric 
factors  occur  in  ratio  pairs.  Consequently,  if  the  volume  of  the  vapor-phase 
sample  bulb  is  chosen  so  that  nu  is  approximately  equal  to  nd.  the  two  sets  of 
manometric  measurements  are  carried  out  under  similar  conditions,  and 
systematic  errors  in  the  numerator  and  denominator  tend  to  cancel.  The 
acceleration  due  to  gravity,  which  is  implicit  in  the  pressures,  need  not  be 
known.  The  volumes  Vs  and  Vb  of  the  bulbs  are  determined  as  a  function  of 
temperature  by  mercury  weighing.  For  the  density  of  mercury  we  use  the 
values  given  by  Bigg,136'  but  again,  as  in  the  pressure  measurements,  any 
possible  systematic  error  would  cancel.  Even  the  effect  of  small  amounts  of 
gas  impurities  tends  to  cancel  if  the  solubilities  of  the  impurities  are  not  too 
different  from  that  of  the  gas  being  studied. 

In  Eq.  (17)  the  uncertainty  in  the  gas  constant  R  is  approximately  30  ppm 
(Sengers  ct  a/.'371),  and  the  errors  in  M  and  p  are  less  than  5  ppm.  The  equilib¬ 
rium  temperature  is  measured  with  a  quartz  crystal  thermometer  (Hewlett- 
Packard  Model  2801  A).'38’  We1391  have  found  that  with  appropriate 
calibration,  it  will  measure  temperature  relative  to  1PTS-68  with  an  absolute 
accuracy  of  probably  0.003  K  and  no  worse  than  0.005rK  (approximately 
1 5  ppm).  In  both  the  large  constant-temperature  baths  designed  for  this  work, 
control  varies  from  +0.001  K  at  lower  temperatures  to  ±0  003  K  (10  ppm) 
at  60  C.  The  crystal  thermometer  probe  is  located  near  the  equilibrator  in 
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the  bath,  and  tests  show  the  maximum  difference  in  temperature  between  the 
bath  and  the  solution  in  the  equilibrator  is  no  greater  than  0.002  K.  (7  ppm). 
For  oxygen  between  0  and  60  C,  an  error  of  no  more  than  80  ppm  is  intro¬ 
duced  by  assuming  there  is  no  difference  between  IPTS-68  and  the  thermo¬ 
dynamic  temperature  scale. 

The  liquid-phase  bulbs  have  volumes  of  800  to  900  cm1  with  an  un¬ 
certainty  less  than  10  ppm.  As  discussed  above,  the  smaller  vapor-phase 
bulbs  are  chosen  in  si2e  according  to  both  the  gas  and  the  temperature  to 
make  nu  approximately  equai  to  na.  For  oxygen  the  maximum  error  in  l'„ 
varies  from  20  ppm  at  0  C  to  40  ppm  at  60  t\  The  precision  of  the  manometric 
determinations  of  naR  and  n&R  depends  upon  the  amount  of  gas  and  there¬ 
fore,  for  approximately  fixed  FR.  upon  both  the  gas  and  the  temperature.  Each 
gas  sample  is  measured  twice,  with  the  gas  expanded  and  recompressed 
between  analyses  to  make  them  independent.  Twice  the  standard  deviation 
of  the  percentage  difference  of  the  two  measurements  from  each  other  for  a 
number  of  samples  provides  a  useful  indication  of  the  manometric  precision 
For  oxygen  it  ranges  from  40  ppm  (0  C)  to  80  ppm  (60  C). 

From  these  estimates  we  might  expect  measurements  of  T.  Fs,  rb, 
nuR.  and  naR  to  lead  to  a  total  random  analytical  error  in  A  of  roughly  0.01  % 
at  low  temperatures  and  0.02^,  at  60  C.  The  total  systematic  analytical  error 
from  R,  M.  />,  and  T.  assuming  the  worst-case  situation,  should  be  less  than 
0.01 2%. 

The  extraction  and  gas  processing  procedures  have  been  tested  exhaus¬ 
tively.  For  example,  although  the  partial  pressures  in  the  drying  manifolds 
do  not  exceed  the  saturated  vapor  pressure  of  most  gases  at  the  temper¬ 
ature  of  liquid  nitrogen,  tests  are  carried  out  with  each  gas  to  show  that 
it  is  not  condensed  in  the  trap.  No  gas  is  mechanically  captured  in  the  dry  ice 
trap  by  freezing  water  vapor.  We  have  found  that  drying  successive  slugs  of 
gas,  which  transfers  relatively  little  water  to  the  trap,  is  preferable  to  pumping 
straight  through  and  thereby  transferring  large  quantities  of  water.  In  another 
test,  measured  quantities  of  gas  in  a  manometer  bulb  are  transferred  to  the 
input  of  the  extraction  vessel,  where  the  liquid-phase  sample  bulb  normally  is 
attached.  With  gas-free  water  in  the  extraction  vessel,  the  gas  is  let  into  the 
vessel  and.  after  stirring,  the  gas  is  extracted  in  the  normal  manner.  When  it  is 
measured,  it  agrees  with  the  initial  amount  In  a  different  test  for  possible 
losses  of  gas  in  the  processing,  a  measured  quantity  of  gas  in  the  manometer 
is  pumped  with  an  extraction  manifold  Toepler  pump  into  its  storage  bulb. 
When  the  gas  is  pumped  back  into  the  manometer  with  the  Toepler  pump  in 
the  manometric  system  and  remeasured,  no  signilicant  change  in  the  number 
of  moles  is  detectable  after  one  such  cycle,  even  if  the  gas  is  left  in  the  storage 
bulb  overnight.  Repealed  cycling  of  this  kind  does  show  that  for  oxygen  an 
average  loss  of  approximately  30  ppm  per  cycle  occurs,  but  even  if  this  were 
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detectable  in  the  manometer,  the  errors  in  the  dissolved  and  undissolved  gas 
sample  determinations  would  at  least  partly  cancel  as  discussed  earlier. 

In  summary,  with  good,  clean,  high-vacuum  technique,  including  simple 
standard  vacuum  tests  before  and  after  processing  a  gas  sample,  the  whole 
operation  of  extracting,  handling,  and  measuring  the  gas  is  straightforward, 
reliable,  and  quantitative. 

The  equilibrator  has  been  tested  thoroughly  in  a  variety  of  ways.  For  the 
large  volume  of  liquid  involved  (approximately  1.5  liters),  the  rate  of  equili¬ 
bration  is  relatively  rapid.  Although  the  design  of  the  filling  arrangement  does 
not  permit  accurate  short-term  kinetic  studies,  equilibrations  over  longer 
periods  of  time  suggest  that  at  5  C  and  the  normal  pumping  speed  of  800 
cm1  min  ',  the  half-time  for  equilibration  is  4  or  5  min.  This  corresponds  to 
approximately  2  to  2.5  cycles  of  the  liquid  through  the  system.  Equilibration 
periods  of  from  4  to  1 1 6  h  yield  identical  results  within  analytical  error, 
indicating  both  achievement  of  equilibrium  and  absence  of  vacuum  problems. 
Mass  spectrometric  tests  show  no  growth  of  impurities  in  the  gas  during 
either  equilibration  or  later  storage  in  glass  sample  tubes  before  the  isotopic 
analyses.  (It  should  be  noted,  however,  that  silicone  vacuum  grease  is 
imperative  when  storing  oxygen  for  even  a  short  time  interval.  Other  greases 
are  oxidized,  with  accompanying  fractionation  of  the  oxygen  isotopes.) 
Equilibration  is  slower  at  high  temperatures  because  the  fourfold  decrease  of 
the  viscosity  .  educes  both  the  pumping  action  of  the  pump  and  mixing  within 
the  liquid-phase  sample  bulb.  An  overnight  equilibration  period  of  20  h  has 
been  found  to  be  convenient  for  all  but  the  highest  temperatures,  where  44  h 
is  used. 

Although  varying  the  pumping  speed  over  the  range  500  to  900  cm3- 
min  1  does  affect  the  rate  of  equilibration,  it  does  not  change  the  results  for 
A.  The  pumping  speed  is  chosen  primarily  to  produce  very  stable,  non- 
turbulent  flow.  In  normal  operation  it  is  so  smooth  that  it  is  difficult  to 
observe  the  flow  of  water  over  the  inner  sphere.  Approach  to  equilibrium 
from  the  supersaturated  state  (achieved  by  first  equilibrating  at  a  lower 
temperature)  leads  to  values  for  A  that  are  identical  with  those  from  the 
undersaturated  state.  Several  equilibrator  units  have  been  used,  including 
ones  without  the  expansion  cylinder  and  one  using  greased  stopcocks  on  the 
liquid-phase  sample  bulb.  Although  all  give  the  same  results  within  analytical 
error,  the  greased  unit  is  not  used  because  the  grease  makes  it  difficult  to 
keep  the  system  clean. 

All  tests  have  shown  the  equilibrator  to  be  extremely  reliable.  The 
design  eliminates  the  thin-film  problems  discussed  in  Appendix  A.  and  yields 
values  of  A  for  the  bulk  liquid.  We  believe  that  there  are  no  systematic  errors 
associated  with  it  and  that  the  accuracy  of  the  results  is  determined  by  the 
precision  of  the  analytical  measurements. 
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2.2.  Isotopic  Fractionation 

The  equilibrium  isotopic  fractionation  factor  a  for  the  solution  of  3402 
and  3202  is  defined  as  the  abundance  ratio  of  the  two  molecules  in  the 
dissolved  gas  to  the  corresponding  ratio  in  the  gas  above  the  solution.  In  the 
notation  employed  above. 


_  34nd/32«d  =  U 
~  34nJ3*na  ru 


(18) 


With  algebraic  rearrangement  and  the  use  of  Eq.  (17).  this  becomes 


a  = 


22k 

3*k 


(19) 


[Negligible  error  is  introduced  by  assuming  that  the  factor  (I  +  APt)  is  the 
same  for  the  two  molecular  species.]  Consequently,  measurement  of  the 
isotopic  fractionation  factor  yields  the  ratio  of  the  Henry  coefficients  for 
the  two  gases. 

Mass  spectrometric  techniques  make  it  possible  to  determine  a  with  very 
high  accuracy.  Ion  beams  produced  from  the  two  isotopic  molecule  are 
simultaneously  collected  in  separate  Faraday  cups,  and  a  quantity  P,  pro¬ 
portional  to  their  ratio,  is  measured  with  a  voltage  divider  in  a  null  arrange¬ 
ment.  Then  by  feeding  the  undissolved  and  dissolved  gases  alternately  into  the 
mass  spectrometer  ion  source,  a  direct  determination  of  the  difference 
±P  =  PA  -  Pu  is  made.  From 

8  =  a  —  1  =  -u-  =  ^  (20) 

ru  ”u 

8  and  a  can  be  calculated.  Systematic  effects  like  ion-source  discrimination 
tend  to  cancel,  and  with  appropriate  techniques  the  precision  and  accuracy 
usually  are  limited  by  noise.  For  the  work  reported  here,  the  values  of  8  lie 
within  the  range  0.5  to  0.9  per  mil  and.  from  repeated  mass  spectrometric 
measurements  on  the  same  pair  of  gas  samples  and  on  gases  from  separate 
equilibrations  at  essentially  the  same  temperature,  the  precision  is  better  than 
0.02  per  mil.  This  means  that  the  ratio  a  =  32kfa*k  can  be  measured  with 
an  error  less  than  20  ppm,  and  the  uncertainty  in  the  values  for  3,A  is  deter¬ 
mined  by  the  accuracy  of  32A'  rather  than  by  the  mass  spectrometry. 

For  the  isotopic  measurements  we  used  a  modified  version  of  the  single- 
focusing.  60  magnetic  sector  mass  spectrometer  described  by  Nier.140’  It  is 
equipped  with  a  dual  system''1"  for  rapid  interchange  of  samples  into  the  ion 
source,  and  the  mass  spectrometer  tube  is  designed  to  permit  null  measure¬ 
ments  for  pairs  of  ions  with  mass  ratios  up  to  3/2  by  insertion  of  the  appro¬ 
priate  dual  collection  system.  To  measure  small  isotopic  effects  of  the  kind 


668 


Benson,  Krause,  and  Paterson 


reported  here,  it  is  important  that  the  mass  spectrometric  analyses  be  made 
twice,  one  immediately  after  the  other,  with  the  gas  samples  exchanged 
between  the  two  sample  systems  between  analyses.  This  compensates  for 
possible  differences  in  the  fractionation  caused  by  the  capillary  leaks,  and  for 
the  minute  differences  in  the  background  contaminations  from  the  sample 
systems.  Although  the  latter  effect  normally  is  small,  we  have  found  that, 
even  with  a  very  "clean”  spectrometer,  it  usually  is  the  more  significant 
source  of  potential  error.  The  average  correction  in  the  oxygen  measure¬ 
ments  to  be  reported  here  was  only  0.007  per  mil,  but  in  isotopic  analyses  of 
some  other  gases,  corrections  up  to  0.05  per  mil  have  been  required. 

3.  RESULTS  AND  DISCUSSION 

3.1 .  Data  j 

According  to  Eq.  (17),  at  fixed  temperature,  k'  should  be  a  linear  function  I 

of  Pl  with  intercept  k  and  slope  kX.  We  have  made  sets  of  measurements  of 
k'  vs.  Py  for  oxygen  at  each  of  four  temperatures,  0.  1 5.  35,  and  55‘C.  Two  of 

these  are  illustrated  in  Fig.  2,  where  it  is  evident  that  within  experimental  ■ 

error  the  relationships  are  linear  over  the  indicated  pressure  ranges.  Although 

the  value  of  A  at  each  temperature  can  be  calculated  directly  from  the  ratio  ' 

of  the  slope  of  the  line  to  its  intercept,  it  is  useful  to  replot  the  data  as  k '/k  vs.  1 

Pi  because  the  slopes  of  the  resulting  lines  are  the  values  of  the  corresponding  , 
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A's.  This  is  shown  in  Fig.  3,  where  it  is  clear  that  A  decreases  as  the  temperature 
increases.  Furthermore,  a  graph  of  A  vs.  t  is  linear  according  to 

A  =  0.002805  -  0.0000375f  (21) 

where  t  is  in  degrees  Celsius.  The  reason  for  this  empirical  result  is  not 
obvious.  Recall  that  if  one  uses  the  Lewis  and  Randall  rule,  A  =  y  +  0t  +  6. 
Furthermore,  •/  is  small.  Therefore,  because  values  for  8  calculated  from  the 
second  virial  coefficient  for  oxygen  between  0  and  60  C,  given  by  Sengers 
cl  al..,3~'  fit  approximately  the  function 

8  =  0.000975  -  1.426  x  10  sr  +  6.436  x  10  et2  (22) 

it  is  surprising  that  A  should  be  linear  in  i.  With  A  known  as  a  function  of  /, 
Eq.  (17)  can  be  used  to  find  A  at  any  temperature  from  a  single  set  of  measure¬ 
ments  at  any  total  pressure.  Experimental  results  for  A  and  S  from  0  to  60  C 
are  given  in  Table  I.  As  noted  earlier,  we  believe  these  values  for  A  range  in 
accuracy  from  approximately  0.01%  at  0  C  to  0.02%  at  60  C,  and  the  values 
for  5  are  probably  good  to  better  than  0.02  per  mil. 

3.2.  The  Variation  of  the  Henry  Coefficient  with  Temperature 

It  is  interesting  to  examine  empirically  how  various  functions  fit  the 
accurate  data  presented  here.  The  results  for  several  types  of  series  are  shown 
in  Table  11.  The  data  have  been  fitted  to  each  of  the  series  by  straightforward 
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Table  O.  Percent  Standard  Deviat.ons  of  k  when  the 
Experimental  Values  for  In  k  and  T  Are  Fitted  to 
Series  of  Various  Forms  and  Numbers  of  Terms 


Form 

Number  of  terms 

3 

4 

5 

6 

l  1  IT 

0.0176 

0.0178 

0.0178 

0.0164 

1  I/In  T 

0.1338 

0.0181 

0.0178 

0.0163 

Valentiner“ 

0.1358 

0.0193 

0.0181 

0.0164 

v  In  T 

0.2167 

0.0197 

0.0180 

0.0163 

1  T 

0.4887 

0.0421 

0.0191 

0.0164 

*  The  Valcntiner  series  is  Eq.  (1). 


matrix-inversion  techniques142’  on  a  time-sharing  terminal  of  a  CDC  Cyber 
74  computer.  In  each  case  we  have  used  from  three  to  six  terms  in  the  ex¬ 
pansion,  and  then  examined  the  resulting  standard  deviations.  It  is  evident 
that  the  power  series  in  1/7"  provides  by  far  the  best  fit.  In  fact,  the  three-term 
expansion  in  1/7"  fully  fits  the  data,  i.e.,  additional  terms  provide  no  improve¬ 
ment  and  the  precision  of  fit  is  comparable  to  the  precision  of  the  measure¬ 
ments.  For  all  the  other  types  of  series  at  least  four  terms  are  required,  and 
further  improvements  are  achieved  by  adding  a  fifth  term.  In  a  related  test 
the  data  have  been  fitted  to  three-term  power  series  in  l/T*,  i.e., 

In  k  =  a0  +  aJT *  +  aJ(T *)2  (23) 

Figure  4  shows  the  variation  of  the  percentage  standard  deviation  of  the 
data  from  each  series  as  a  function  of  the  parameter  .v.  The  minimum  devia¬ 
tion  occurs  with  v  very  nearly  equal  to  unity  (.vmln  =  0.99)! 

The  Valentiner  series,  Eq.  (1),  is  an  integrated  van’t  Hoff  equation  that 
in  one  form  or  another  has  been  employed  for  many  years  to  express  the 
temperature  dependence  of  equilibrium  constants  In  1966,  Clarke  and 
G lew1 231  presented  a  generalized  treatment  of  the  integration  with  an 
excellent  discussion  of  the  evaluation  of  the  errors  in  the  standard  thermo¬ 
dynamic-function  changes  when  Eq.  (1)  is  used  to  fit  experimentally  deter¬ 
mined  equilibrium  constants.  In  their  formalism,  Eq.  (1)  follows  from  an 
expansion  of  the  standard  enthalpy  change  in  a  Taylor  series  in  T  about 
a  reference  temperature  6.  As  such  it  is  undeniably  an  appropriate  expression 
for  fitting  experimental  data.  More  is  claimed  for  this  expression,  however. 
Clarke  and  Glew  state  that  “only  those  equations  equivalent  to  (Eq.  (1)] 
furnish  the  unique,  best  unbiased  estimates  of  the  standard  thermodynamic- 
function  changes  for  reactions."  This  is  too  sweeping.  Thermodynamic  argu¬ 
ments  alone  cannot  lead  one  to  predict  the  otherwise  unknown  temperature 
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Fig.  4.  Variation  of  the  percentage  standard  deviation  of  the  data  from  the  equation 
In  k  =  a0  +  aJT*  +  a2!(T*)2  as  a  function  of  the  parameter  x.  Minimum  deviation 
occurs  at  jr  =  0.99. 

dependence  of  equilibrium  constants.  The  danger  in  prescribing  one  functional 
form  is  that  it  may  discourage  the  use  of  other,  simpler,  and  equally  valid 
forms  which  may  give  more  insight  into  the  underlying  phenomena. 

In  Appendix  B  we  show  that  Eq.  (2).  the  power  series  in  \/T,  may  be 
derived  in  a  manner  that  is  identical  to  that  given  by  Clarke  and  Glew  for 
Eq.  (I),  except  that  a  different  measure  of  temperature,  1  IT,  is  used  in  place 
of  T.  (One  could  argue  that  1/T  is  fundamentally  just  as  good  a  measure  of 
temperature  as  T ,  or  even  better  because  it  occurs  more  naturally  in  statistical 
mechanics.)  Expansion  of  AH'  in  a  Taylor  series  in  \/T  about  the  reciprocal 
reference  temperature  1/0,  use  of  the  thermodynamic  relationships,  and 
algebraic  manipulation  lead  directly  to  Eq.  (2).  The  coefficients  in  this  series 
are  as  simply  related  to  the  standard  thermodynamic  difference  functions 
as  are  the  coefficients  in  Eq.  ( 1 ),  and  the  method  for  determining  the  functions 
and  their  standard  errors  is  completely  analogous.  Thus,  the  type  of  argument 
that  leads  to  Eq.  (I)  does  not  produce  it  uniquely.  Equation  (2)  and  other 
similarly  consistent  possibilities  are  on  equal  thermodynamic  footing  with 
Eq.  (1). 

On  the  basis  of  thermodynamic  consistency,  ability  to  fit  accurate  data, 
and  simplicity,  Eq.  (2),  the  power  series  in  1/T,  is  clearly  the  best  family  of 
functions  proposed  to  date  for  representing  the  variation  of  k  with  T.  Over 
the  range  of  temperature  studied  here,  its  three-term  member. 

In  k  =  a0  +  aJT  +  aJT2 


(24) 
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is  sufficient  to  fit  the  solubility  data  that  presently  are  available.  Values  for  the 
three  coefficients  for  oxygen  are  given  in  Table  VI.  It  remains  to  be  seen 
whether  there  is  a  fundamental  significance  to  Eq.  (24).  The  idea  seems  to  be 
supported  by  the  relationships  between  the  molecular  properties  of  the  noble 
gases  and  the  parameters  of  their  solutions,  which  emerged  from  the  analysis 
of  the  preliminary  data.'281  Below,  we  shall  see  that  the  present  results  for 
oxygen  are  in  accord  with  that  analysis. 

Although  the  solubility  measurements  reported  here  apparently  are  the 
most  accurate  determinations  ever  made  of  the  equilibrium  constant  for  any 
reaction  over  a  wide  range  of  temperature,  their  excellent  fit  to  Eq.  (2)  does 
not  necessarily  imply  that  the  equation  should  be  applicable  to  other  kinds 
of  reactions  where  the  interaction  mechanisms  are  different.  We  have  ex¬ 
plored  the  literature  for  the  best  data  on  equilibrium  constants  for  other  types 
of  systems  such  as  weak  acids  and  bases  in  water.  Most  of  the  data  are  either 
too  imprecise  or  too  limited  in  temperature  range  to  distinguish  between  the 
ability  of  Eqs.  (1)  and  (2)  to  fit  them.  For  the  few  relatively  accurate  results 
available,  viz.  the  work  of  Feates  and  Ives'43’  and  Ives  and  Marsden1441  on 
weak  acids,  Eq  (2)  gives  the  better  fit.  and  again  only  the  three-term  equation 
(24)  is  required,  while  use  of  Eq.  ( 1 )  requires  four  terms.  These  results  suggest 
that  Eq.  (24)  [or.  more  generally,  Eq.  (2)]  may  have  broader  applications 
than  for  gas  solubilities,  but  precise  data  for  a  number  of  other  systems  are 
needed  to  test  this  possibility. 


3.3.  Comparison  with  Previous  Measurements 

The  fifth  column  of  Table  1  gives  the  percentage  deviations  of  the  new 
data  for  the  Henry  coefficient  of  oxygen  from  their  fit  to  Eq.  (24),  the  three- 
term  power  series  in  I 'T.  These  are  illustrated  graphically  in  Fig.  5.  where  the 
new  data  are  represented  by  the  solid  circles  distributed  along  the  reference 
line.  For  comparison  with  probably  the  best  previous  determinations,  we  have 
recalculated  the  experimental  values  of  Montgomery.  Thom,  and  Cock- 
burn,'1'7'  Murray  and  Riley,'22' Klots  and  Benson,"61  and  Carpenter120’  to  give 
values  for  k  with  appropriate  corrections  to  make  them  correspond  as  nearly 
as  possible  to  our  values.  The  results  are  included  in  Fig.  5.  The  very  high 
precision  of  the  new-  results  is  evident,  but  another  observation  is  equally 
significant.  Visual  inspection  indicates  that  the  five  sets  of  results  are  in 
excellent  agreement,  with  the  four  older  sets  of  data  points  scattered  more  or 
less  equally  on  either  side  of  the  smoothed  curve  for  the  new  data.  Table  111 
illustrates  this  in  another  way.  The  rms  percentage  deviation  for  the  new 
data  is  more  than  an  order  of  magnitude  smaller  than  that  for  any  of  the  other 
four  sets  of  data.  The  percentage  deviations  of  the  Klots  and  Benson  set.  the 
Carpenter  set,  and  the  Montgomery  et  al.  set  average  almost  to  zero,  indi- 
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Fig.  5.  Variation  with  temperature  of  the  percentage  deviation  of  the  new  measurements 
for  A,  and  of  the  best  previous  data  for  oxygen,  from  the  equation  In  A  =  3.71814  + 
5595  17  T  -  1049068 T2.  The  smooth  curve  of  Wilhelm  el  al.  (ref.  3)  is  also  shown. 


eating  each  is  equally  scattered  about  the  new  fit,  while  the  values  of  k  from 
Murray  and  Riley  average  0.18%  higher. 

This  agreement  is  very  important,  because  it  suggests  that  the  results  are 
probably  free  from  systematic  errors  well  within  the  random  errors  in  the 
older  data.  On  the  basis  of  our  own  tests,  as  discussed  in  Section  2.1,  we 
believe  the  accuracy  of  the  new  technique  is  determined  by  the  random 
errors,  i.e.,  approximately  0.01  to  0.02%.  but  the  agreement  with  the  best 


Table  UI.  Comparison  of  the  Deviations  of  Data  from  Various  Sources 
with  the  Present  Results 


Source 

V  %  deviation 

n 

(T  deviation)2 j 1,3 

Murray  and  Riley  " 

+  0  180 

0.264 

Klots  and  Benson" 

+  0.042 

0.212 

Carpcntcrr 

-0.024 

0.220 

Montgomery,  Thom,  and  Cockburn* 

-  0.079 

0.215 

Present 

0 

0.017 

°  See  ref.  22.  6  Sec  ref  16. 

'  See  ref.  20.  ■‘See  ref.  17. 
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previous  solubility  measurements  on  any  gas  adds  confidence  to  that  judg¬ 
ment.  This  independent  confirmation  of  the  method  is  significant  also  for 
results  we  shall  report  later  on  other  gases,  because  for  them  the  previous 
data  are  of  lower  quality  than  for  oxygen  and  no  significant  comparison  with 
the  new  results  will  be  possible. 

The  curve  obtained  by  Wilhelm,  Battino,  and  Wilcock<3)  by  fitting 
selected  data  to  a  four-term  series  of  the  Valentiner  type,  Eq.  ( I ),  is  included 
in  Fig.  5.  The  comparison  is  not  significant  with  respect  to  the  merits  of 
Eq.  (1)  relative  to  Eq.  (2)  because  other  less  accurate  data111’  were  used  by 
Wilhelm  et  al.  for  the  high-temperature  region,  but  it  is  clear  that  above  40°C 
their  smoothed  values  become  quite  incorrect. 

3.4.  Thermodynamic  Properties  of  the  Solution 

For  systems  of  this  kind,  the  standard  state  of  the  gas  in  the  vapor  phase 
at  any  temperature  is  usually  defined  to  be  that  in  which  the  gas  has  a  fugacity 
of  1  atm,  while  the  standard  state  for  the  dissolved  gas  at  any  temperature 
is  the  hypothetical  state  found  by  extrapolating  the  line  representing  Henry’s 
law  on  an  /  vs.  .v  graph  to  unit  mole  fraction  for  the  dissolved  gas.  With  these 
choices,  and  Eq.  (24),  the  changes  in  the  partial  molal  thermodynamic 
functions  are 


AC°  =  R[a0T  +  al  +  aJT ] 

(25) 

Atf°  =  /?[<*!  +  2o2/T] 

(26) 

AS0  =  /?[-fl0  +  a2/T2) 

(27) 

A  C;  =  -2  RaJT2 

(28) 

[Note  that  these  differ  in  sign  from  those  given  earlier'28'  because  In  k, 
rather  than  ln(l/Ar),  has  been  used  in  Eq.  (24).]  Values  for  AG°,  A H°,  A S°, 
and  A C°  are  given  at  four  temperatures  in  Table  IV,  with  numbers  calculated 
by  Wilhelm,  Battino,  and  Wilcock'3’  and  the  theoretical  predictions  of 
Pierotti.'45'  The  percentage  standard  deviation  in  our  values  for  In  k  (and, 
therefore,  A G°)  is  0.001 67>-  This  leads  to  a  random  error  in  AC£  equal  to 
0.1 570.  with  a  9570  confidence  level.  As  would  be  expected  from  Fig.  5,  the 
improvements  over  the  Wilhelm  et  al. t3)  results  are  small  except  for  AC*  and 
at  higher  temperatures.  The  Pierotti  values  are  systematically  below  the 
experimental  ones,  with  AC£  smaller  by  approximately  30%  at  298°K. 
Reference  to  Pierotti’s  paper  also  shows  that  the  predicted  rate  of  change  of 
AC°  with  temperature  is  only  1/3  the  experimental  value.  It  should  be  noted, 
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Table  IV.  Comparison  of  Vr  lues  for  A(7\  AW ",  AS°,  and  AC*  with  Litera¬ 
ture  Values1' 


Source 

AG 

(cal-mole  1 

AH 

)  (cal-mole  1 

AS” 

')  (cal-^K  * 1  mole 

AC° 

‘Mcal-'K  ‘-mole'1) 

Present 

5846.0 

283.15  K 
-3613 

—  33.41 

52  03 

Wilhelm  el  at  * 

5847 

-  3594 

-33.34 

47.19 

Present 

6327.4 

298.15  K 
-2872 

-30.85 

46.93 

Wilhelm  et  al.“ 

6328 

-2882 

-30.89 

47.75 

Pierotti,  theoretical” 

6270 

-2423 

-28.5 

32.2 

Present 

6773.3 

313.15  K 
-2201 

-28.66 

42.54 

Wilhelm  et  ai* 

6774 

-2162 

-28.53 

48.31 

Present 

7188.7 

328.I5°K 

-1592 

-26.76 

38.74 

Wilhelm  et  al.a 

7184 

-1433 

-26.26 

48.88 

0  See  ref.  3. 

"  See  ref.  45. 

f  R  =  1.987(7  ca!-°K  ‘-mole1. 


especially,  that  Eq.  (28)  states  that  AC°  varies  inversely  with  the  square  of 
the  absolute  temperature. 


3.5.  Solution  Parameters  and  Molecular  Properties  of  the  Gases 

From  analyses  of  preliminary  experimental  data  for  seven  gases,  Benson 
and  Krause'28’  found  that  when  the  variation  of  k  with  T  was  expressed  by 
Eq.  (24).  relationships  emerged  among  the  solution  parameters  and  the 
thermodynamic  and  molecular  properties  of  the  gases.  This  became  partic¬ 
ularly  apparent  when  Eq.  (24)  was  rewritten  in  the  form 

In  A-  =  -  l)  +  Aa(T±  -  l)2  (24a) 

and  the  data  were  smoothed  by  requiring  that  the  constant  A2  be  the  same  for 
all  the  gases.  (Some  other  parameters  are  defined  in  terms  of  a0.  0i,  and  a2 
in  Table  V.)  Benson  and  Krause  found  that  “although  the  A2  smoothing 
slightly  increases  some  values  and  decreases  others,  for  ail  parameters  the 
patterns  of  the  numbers  are  skewed  very  little  in  any  systematic  way." 
Furthermore,  each  of  the  parameters,  a0,  «i.  o2,  Tu  A3,  and  others,  became  an 
essentially  linear  function  of  (*/As)!  2,  the  square  root  of  the  force  constant 
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Table  V.  Definitions  of  Some  Parameters 


7)  =  -  -taoUj)1  2)  A o  =  do 

As  =  -(a?  -  4a0«a)’  /),  =  a,ITt 

M  —  [-a?  +  4u0aa]/4d2  A2  =  a2iTl 

7'm  =  —  2d2/di  /40  +  rli  4-  ,4Z  =  0 


-  7)  is  Ihe  absolute  temperature  at  which  the  Henry  coeflicient  A 
hypothetically  would  be  unity,  if  the  systems  hypothetically  had 
the  same  properties  outside  the  experimental  temperature 
range. 

b  47  ~  In  Amux- 

c  7\»  is  the  temperature  for  maximum  A. 

of  the  gas.  They  concluded,  "The  implication  is  strong  that  with  more  accurate 
measurements  A2  and  (Afp),  would  be  found  to  be  the  same  for  all  the  gases." 

A  partial  test  of  this  suggestion  is  possible  from  a  comparison  of  the 
values  for  the  parameters  obtained  here  with  the  corresponding  values  pre¬ 
dicted  from  the  preliminary  data  smoothed  on  the  assumption  that  A2 
should  be  common  to  all  the  gases.  These  are  shown  in  Table  VI,  together 
with  the  preliminary  '  raw”  values.  The  improvement  achieved  by  the  A2 
smoothing  process  is  striking,  and  it  lends  support  to  the  idea  that  Eq.  (24a) 
is  a  universal  solubility  equation  in  the  sense  that  A2  is  common  to  all  the 
gases  with  only  7",  and  A3  being  characteristic  of  the  specific  gas. 


Table  VI.  Comparison  of  the  Values  for  Various  Parameters  with  the 
Smoothed  and  Raw  Values  Obtained  in  the  Preliminary  Work  (Reference  28) 


a i  ~<?2 

Source  a0  =  A0  (°K  x  10  3)  ('K!  x  10  5|  7\  {“ K) 


Present 

3.71814- 

5.59617- 

10.49668“ 

168.667" 

Preliminary,  smoothed 

3.7671 

5.5869 

10.5075 

168.85 

Preliminary,  raw 

4.0605 

5.4167 

10.261 

168.22 

Source 

.4, 

~  ^2 

~  -4  3 

M 

T„  (°K) 

Present 

33.179 

36.897 

11.177 

375.14 

Preliminary,  smoothed 

33.088 

36.855 

40.622 

11.193 

376.15 

Preliminary,  raw 

32.200 

36.261 

40.321 

1 1 .209 

378.87 

‘  The  three  present  primary  parameters,  a0,  Oi.  and  a2.  and  7\  individually  are  not 
accurate  to  the  number  of  figures  given,  but  they  are  coupled  by  the  fitting  process  and 
it  is  useful  to  carry  an  extra  figure  for  internal  consistency  in  parameter  conversion. 
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Equation  (24a)  expresses  k  in  terms  of  the  “scaled  temperature"  TjT ,. 
The  equation  also  may  be  written  as 

In  A-  =  A„  +  +  ^(y)2  (24b) 

which  is  similar  in  form  to  Eq.  (24),  but  here  the  third  term  on  the  right,  and 
therefore  ACP  =  -  2RA.,(TJT )2,  would  be  the  same  for  all  gases  at  the  same 
scaled  temperature.  (Note  that  because  T  =  7\  when  k  =  I,  A0  +  Ax  +  A2 
=  0,  and  Eq.  (24b)  involves  only  three  independent  constants.) 

3.6.  Values  for  k  at  Temperatures  above  60°C 

Values  for  k  at  high  temperatures,  well  beyond  the  range  of  our  deter¬ 
minations,  are  needed  for  many  purposes.  Although  measurements  that  have 
been  reported  for  temperatures  above  100  C  are  much  less  accurate,  we  have 
attempted  to  relate  them  to  our  work.  The  result  is  only  partly  successful, 
but  at  least  for  the  present  it  provides  both  the  best  and  easiest  way  to  calculate 
values  for  k  at  elevated  temperatures. 

To  test  the  ability  of  our  data  and  Eq.  (24)  to  give  correct  values  when 
extrapolated  to  higher  temperatures,  the  data  for  four  different  temperature 
ranges,  0-20,  0-30.  0-40,  and  0-60  C.  have  been  fitted  separately  to  Eq.  (24) 
and  then  used  to  predict  k  at  temperatures  above  the  range  of  each  set  of  data. 
The  results  are  shown  in  the  first  four  rows  of  Table  VII.  (The  "extrapolated 
value  ”  at  60  C  for  the  0-60'C  range  is  underlined  to  indicate  it  is  not  extrap¬ 
olated.)  The  predicted  or  extrapolated  values  from  the  data  for  0-20.  0-30, 
and  0  40  C  agree  extremely  well  with  the  measured  value  at  60  C.  For 
example,  the  0-30  C  results,  when  extrapolated  to  twice  the  temperature 
range  of  the  data,  disagree  with  the  measured  value  by  only  0.06%. 

The  most  extensive  measurements  of  oxygen  solubility  at  high  tempera¬ 
tures  and  pressures  have  been  those  carried  out  at  the  Battelle  Memorial 
Institute  and  included  in  the  report  by  Stephan.  Hatfield,  Peoples,  and 
Pray.'1"’  They  defined  k  as  the  ratio  of  the  partial  pressure  of  the  gas  in 
pounds  per  square  inch  absolute  to  the  solubility  in  milliliters  of  gas  at  STP 
per  gram  of  solution,  in  the  limit  as  the  partial  pressure  approaches  zero. 
In  the  fifth  row  of  Table  VII  we  have  reexpressed  their  results  in  atmospheres, 
but  their  experimental  method  does  not  permit  calculation  of  a  k  that  w  ould 
correspond  exactly  to  the  Henry  coefficient  defined  in  Eq.  (3).  Despite  the 
lack  of  equivalence  of  the  two  A's.  it  is  useful  to  compare  the  Battelle  experi¬ 
mental  results  with  the  values  obtained  by  extrapolating  the  Eq.  (24)  fit  to 
our  complete  set  of  data  for  0-60  C.  At  100  C  the  numbers  differ  by  1.9%, 
which  is  within  the  uncertainty  of  their  experimental  numbers.  At  higher 
temperatures,  however,  the  difference  becomes  progressively  greater. 
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The  significance  of  these  differences  is  not  clear,  hut  it  i*  interesting  to 
note  that  application  of  Eq  (2)  with  four  terms  to  a  combination  of  the  data 
of  Stephan  ct  al.  for  100-288  C  with  our  data  for  0  60  C  ( with  equal  weight¬ 
ing  for  all  points  despite  the  disparity  in  precisions)  tits  all  the  Stephan  ct  al. 
data  with  a  rms  deviation  of  2  0%.  which  is  within  their  experimental  error 
Refer  to  the  sixth  row  in  Table  VII  Between  0  and  60  C  the  deviations  con¬ 
siderably  exceed  the  uncertainties  in  our  measurements,  but  they  are  never 
greater  than  0.7"~o.  The  fit  to  the  high-temperature  data  is  better  than  that 
given  by  the  more  complex  Himmelblau'25''  equation,  and  calculation  of  A 
for  a  given  T  is  much  simpler. 

Until  more  accurate  data  become  available  for  temperatures  above  60  C. 
we  recommend  on  the  basis  of  the  results  above  that  v  alues  for  A  be  calculated 
from  one  of  two  equations,  depending  upon  the  temperature  range  involved: 

(a)  Tor  very  accurate  values  in  the  temperature  range  0  100  C: 

Ink  =  2.71814  +  5596.17/7'  -  1049668,-T2  (29) 

where/ =  A. v  (3) 

(b)  For  values  above  100  C': 

In  A  =  -4.1741  +  1.3104  x  1077"  -  3.4170  x  l06/r2 

+  2.4749  x  l08/r'  (30) 

where  p  =  A.v,  and  p  is  a  partial  pressure  (in  atmospheres)  of  the  gas  in  the 
vapor  phase,  calculated  by  subtracting  the  water  vapor  pressure  given  in  the 
steam  tables  from  the  total  pressure,  and  x  is  the  mole  fraction  of  gas  in 
the  solution.  Within  the  experimental  errors,  this  corresponds  approximately 
to  the  Stephan  ct  al.  results,  on  which  Eq.  (30)  is  based. 

3.7.  Isotopic  Fractionation  during  Solution 

In  column  6  of  Table  I  the  measured  values  for  8  range  from  approxi¬ 
mately  0.85  per  mil  («  =  1.00085)  at  low  temperatures  to  0.55  per  mil 
(«  =  1.00055)  at  60  C.  In  Fig.  6.  these  results  are  plotted  as  In  a(  x  I03)  vs. 
1/7".  together  with  the  data  from  Klots  and  Benson'10’  and  Kroopnick  and 
Craig  '"  Note  that  the  ordinate  in  Fig.  6  is  essentially  equivalent  to  8 
(expressed  in  per  mil),  because 

In  «  =  ln(l  +  8)  =  8  -  82/2  +  83/3  +  •  •  • 

Consequently,  since  8  <  0.001  for  oxygen,  82/2  and  higher-order  terms  are 
negligible,  and  to  a  very  good  approximation 
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Fig.  6.  Variation  of  the  logarithm  of  the  fractionation  factor  a  for  oxygen  with  the 
reciprocal  of  the  absolute  temperature.  Results  from  previous  work  arc  also  shown. 


The  high  precision  and  relatively  wide  temperature  range  of  the  new 
measurements  make  possible  an  attempt  to  explore  the  variation  of  a  with 
temperature.  From  Eqs.  (19)  and  (24), 

In  a  =  In  3ik  -  In  34k 

=  ( 33a0  -  34<z0)  +  (32*t  -  34«i )/T  +  C2a2  -  3'a2)/T2 
or 

In  «  =  Aa0  +  Aa,/F  +  Ao2/r2  (31) 

which  is  similar  in  form  to  functions  considered  by  Bigeleisen,'47'  Bardo  and 
Wolfsberg,1481  Rolston  el  a/., 1481  van  Hook  and  Phillips.1501  and  others,  in 
discussing  isotope  effects.  Table  VIII  summarizes  the  results  of  fitting  the 
data  to  various  combinations  of  terms  in  Eq.  (31).  The  standard  deviation  of 
In  a  for  all  forms  is  0.017  per  mil,  which  is  comparable  to  the  estimated  pre¬ 
cision  of  the  measurements.  Clearly,  on  the  basis  of  these  oxygen  data  alone, 
it  is  not  possible  to  choose  one  of  the  functions  as  preferable.  (Although  the 
standard  deviations  are  the  same,  the  deviations  at  any  specific  temperature 
may  differ  slightly.  The  fit  used  in  obtaining  column  7  of  Table  I  is  the  three- 
term  fit.  No.  I  in  Table  VIII.) 

If  desired,  3*k  at  any  temperature  may  be  calculated  from 

In  3*k  =  3.71937  +  S59S.45/T  -  104962 4/T2  (32) 
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Table  VIII.  Fits  of  the  Isotopic  Fraction¬ 
ation  Factor  a  to  Functions  of  the  Form 
In  a  =  A a0  +  A a,/T  +  A aJT2,  for  Various 
Combinations  of  Terms 


(1)  In  a  =  —0.0012263  +  0.72466, T  -  44.398/T2 

(2)  In  a  =  -0.00072951  +  0.42696/J" 

(3)  In  o  =  -0.00001621  +  63. 609/ T2 

(4)  In  a  =  -  0.0092203/7'  +  64.9I4/7'2 

(5)  In  a  =  62.223/7’“ 


“  In  all  cases  the  standard  deviation  of  both  In  a  and 
6  is  0.017V 


where  the  coefficients  have  been  obtained  from  the  relationships  34o,  - 
32a,  -  A  a..  (Again,  the  three-term  function  for  In  a  has  been  used.) 

If  it  is  assumed  that  the  isotopic  fractionation  is  due  to  a  difference  in 
the  zero-point  energies  of  the  two  types  of  oxygen  molecules  in  solution,  then 
a  calculation  like  that  of  Muccitelii  and  Wen  <34’  leads  to  a  value  of  approxi¬ 
mately  2.5  A  for  the  size  of  the  solvent  cage.  For  the  solution  of  hydrogen, 
Muccitelii  and  Wen  obtained  approximately  3  A.  It  might  be  argued  that  a 
smaller  value  should  be  expected  for  oxygen  than  for  hydrogen,  because  the 
larger  oxygen  molecule  would  make  the  effective  size  of  the  cavity  smaller, 
but  additional  precise  measurements  on  other  gases,  and  further  theoretical 
studies  of  quantum  effects  in  aqueous  solutions,  are  needed  before  our 
understanding  becomes  more  than  qualitative. 

It  is  evident  that  differences  between  3202  and  3402  for  the  partial  molal 
function  changes  may  be  calculated  very  simply  from  Eqs.  (25)  to  (28)  by 
substituting  for  a0.  a,,  and  a2  the  coefficient  differences  A  an.  A  a,,  and  Aa2, 
respectively.  For  oxygen,  of  course,  the  difference  in  the  change  in  heat 
capacity  cannot  be  determined  from  the  present  data,  because  reference  to 
the  values  for  A«2  in  Table  Vill  shows  hat  the  difference  in  AC,  could  be 
zero,  negative,  or  positive,  depending  upon  the  function  chosen  for  In  a. 
On  the  other  hand,  the  values  shown  in  Table  IX  for  the  other  differences  in 
thermodynamic  changes  are  independent  of  the  function  chosen.  The  values 
for  (A/7  )/A( A/7  )  and  (A.5  )/A(AS  )  are  interesting  because  they  indicate 
that  substitution  of  an  3,02  molecule  for  an  3202  molecule  in  solution  in¬ 
volves  a  change  in  enthalpy  with  relatively  little  change  in  entropy. 

4.  CONCLUSION 

A  new  method  has  been  presented  for  determination  of  the  Henry 
coefficient  k  and  the  isotopic  fractionation  factor  a  of  gases  dissolved  in 
liquids.  For  oxygen  the  precision  and  probable  accuracy  of  the  measurements 
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Table  IX  Differences  at  298.1 5°K  between 
3202  and  3402  for  the  Partial  Molal  Thermo¬ 
dynamic  Function  Changes  between  the 
Usual  Standard  States 


A(A<7  )  =  0.418  cal-mole'  \G°IM±G')  =  15,100 

A(A/7  )  =  0.85  cal-mole  “ '  A/7/A(A/7  )  =  3,400 

A(  AS  )  =  0  0014  cal-  K.  '-mole  1  ASTAfAS  )  =  21,400 


‘A(AC  )  =  (A(7  ).„  -  (AC  etc. 

of  A  are  an  order  of  magnitude  better  than  previous  results  on  any  gas  This 
very  high  precision  has  shown  that  the  best  form  for  the  variation  of  A  with 
r,  and  perhaps  for  equilibrium  constants  of  other  systems,  is  the  three-term 
power  series  in  IT  proposed  by  Benson  and  Krause  Their  suggestion  that 
the  equation  may  be  rewritten  as  a  universal  solubility  equation  is  supported 
by  the  new  oxygen  results.  The  isotopic  effects  are  small,  but  very  precise 
analyses  have  provided  the  first  quantitative  determination  of  the  variation 
of  n  with  temperature  for  oxygen.  Accurate  values  have  been  presented  for  the 
partial  molal  thermodynamic  function  changes  between  the  liquid  and  vapor 
phases,  including  A Cv  which  varies  inversely  with  the  square  of  the  absolute 
temperature.  The  isotopic  measurements  on  oxygen  have  shown  that  sub¬ 
stitution  of  an  3402  molecule  for  an  3202  molecule  in  solution  involves  a 
change  in  enthalpy  with  a  relatively  small  change  in  entropy.  They  are 
qualitatively  consistent  with  the  idea  that  the  fractionation  is  the  result  of  a 
difference  of  zero-point  energy  of  the  two  molecules  in  solution  A  new 
procedure  has  been  suggested  for  calculating  values  for  A  at  high  temperatures 
until  more  precise  measurements  become  available. 
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APPENDIX  A 

Although  the  equilibrator  employed  in  the  preliminary  measurements 
reported  by  Benson  and  Krause  1281  is  no  longer  used,  a  brief  description  of  the 
experience  with  it  should  be  reported  because  it  suggests  the  possibility  that 
solution  of  a  gas  in  a  very  thin  film  of  liquid  may  be  influenced  by  the  pres¬ 
ence  of  the  surface  supporting  the  film. 

A  drawing  of  the  unit  is  shown  in  Fig.  A-l.  It  consisted  essentially  of  an 
upper  chamber  A,  a  lower  chamber  B,  and  a  connecting  tube  C.  The  total 
volume  of  the  system  was  constant,  but  its  distribution  between  A  and  B 
depended  upon  the  position  of  the  O-ring  sealed  piston.  The  piston  and 


?  f? 

I  | 


Fig.  A-f.  Drawing  of  the  equilibrator  used  by  Benson  and  Krause  (ref.  28)  for  the 
preliminary  measurements. 
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cylinder  were  made  of  Type  316  stainless  steel.  Two  kinds  of  “heads”  for 
chamber  A  were  used.  One  was  blown  from  4-in. -diameter  “Pyrex”  pipe  and 
included  two  6  mm  stopcocks  at  the  top.  The  other  was  constructed  of  Type 
316  stainless  steel  with  specially  designed  stainless  steel  bellows  valves.  The 
O-rings  were  made  of  “  Viton-A,"  and  those  in  contact  with  the  solution  were 
ungreased. 

With  the  piston  in  approximately  the  position  shown,  chamber  A  was 
filled  with  degassed  water  and  then  closed  off.  A  vapor-phase  sample  bulb, 
with  a  stopcock  and  ground  joint  at  each  end,  was  attached  to  the  top  of  the 
equilibrator,  and  the  desired  gas  put  into  the  bulb,  tube  C,  and  chamber  B. 
The  system  was  mounted  vertically  and  fully  immersed  in  a  constant- 
temperature  bath,  the  chambers  were  interconnected,  and  the  piston  drawn 
down  almost  to  the  bottom.  The  support  frame  for  the  equilibrator  was 
tilted  to  approximately  45°.  and  the  unit  was  set  into  rotation  about  its 
cylindrical  axis  to  bring  a  continuously  changing  film  of  liquid  along  the  wall 
of  A  into  contact  with  the  gas  in  essentially  the  same  manner  as  that  described 
by  Green,151’  who  attributed  the  idea  to  Jacobsen.1521  After  the  desired  interval 
of  time,  the  vapor-phase  sample  bulb  was  closed  to  isolate  a  known  volume 
of  gas  at  the  temperature  of  the  bath.  Then,  in  stages,  the  piston  was  raised 
and  the  tilt  reduced,  with  rotation  continued  between  adjustments,  until  the 
unit  was  upright  and  the  lower  flange  of  the  piston  was  firmly  against  the 
shoulder  stop  on  the  cylinder.  During  this  process  the  liquid  level  rose  and 
displaced  gas  from  A  to  B  without  trapping  bubbles  at  the  wall-liquid 
boundary  and  without  pressure  change  within  the  system.  With  A  closed  off 
at  the  top,  an  accurately  known  volume  of  the  solution  was  isolated.  Sub¬ 
sequent  gas  extraction  and  drying  procedures  were  similar  to  those  described 
in  Section  2.1. 

After  eliminating  some  initial  problems  with  the  groove  design  for  the 
upper  O-ring  on  the  piston,  the  equilibrator  gave  very  precise  results  which 
permitted  detection  of  relatively  small,  but  'ignificant,  systematic  effects 
apparently  associated  with  the  thin  film.  Summarized  briefly,  the  following 
characteristics  were  observed: 

(1)  At  all  the  rotation  speeds  that  were  employed,  the  kinetics  of 
equilibration  were  determined  by  the  rate  of  formation  of  film  on  the 
wall  of  A. 

(2)  In  the  early  measurements  the  kinetics  for  a  given  speed  of  rotation 
seemed  to  follow  the  Adeney-Becker153'  growth  curve  (attributed  to 
Carlson (54)  by  Montgomery  el  al.al').  Later,  however,  semilog  plots 
for  long  equilibration  periods  revealed  a  dual  half-time.  Typically, 
for  example,  at  low  temperatures  an  initial  half-time  of  15  min 
shifted  to  a  half-time  of  5  h  when  the  water  became  nearly  saturated. 
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This  seemed  to  imply  a  different  control  mechanism  near  the  end  of 
the  process. 

(3)  As  the  temperature  was  increased,  the  rate  of  equilibration  decreased 
markedly.  This  is  consistent  with  the  very  large  decrease  in  the 
viscosity  of  water  at  high  temperatures,  which  would  lead  to  much 
thinner  films  on  the  surface  and  reduce  mixing  in  the  bulk  of  the 
liquid. 

(4)  Approach  to  equilibrium  from  both  the  undersaturated  and  super¬ 
saturated  states  showed  that,  near  equilibrium,  the  same  concentration 
was  reached  from  both  initial  states,  but  that  the  common  value  for 
the  concentration  increased  and  approached  the  equilibrium  value 
with  the  longer  half-time.  It  is  clear  that  this  traditional  test  for 
equilibrium  is  not  conclusive  unless  the  pairs  of  measurements  are 
shown  to  give  the  same  results  for  several  different  equilibration 
periods.  Later,  after  results  were  available  from  the  equilibrator 
described  in  Section  2.1,  it  became  evident  that,  in  the  very  long 
equilibration  periods,  the  thin-film  equilibrator  could  yield  solutions 
apparently  slightly  supersaturated  relative  to  the  bulk  liquid. 

(5)  The  value  of  k  depended  upon  whether  the  liquid  film  was  in 
contact  with  a  stainless  steel  wall  or  a  glass  wall  in  chamber  A.  For 
example,  with  both  oxygen  and  argon  at  2 °C.  k  was  greater  by 
approximately  O.I70  with  stainless  steel. 

(6)  Tests  showed  that  these  effects  were  not  caused  by  contamination 
from  either  real  or  virtual  vacuum  leaks,  by  incomplete  mixing  in 
chamber  A,  or  by  vacuum  grease  in  the  stopcock  at  the  top  of 
chamber  A  in  the  glass  head. 

After  the  equilibrator  in  Fig.  1  proved  to  give  excellent  results  without 
any  of  the  above  problems,  no  further  attempts  were  made  to  pursue  the 
source  of  the  difficulties.  The  observations  are  consistent,  however,  with  the 
hypothesis  that  the  properties  of  a  thin  film  clinging  to  a  wall  of  metal  or 
glass  may  be  different  from  those  of  the  bulk  liquid.  Adsorption  of  gas  on  the 
large  surface  area  of  chamber  A  may  also  have  played  a  role,  because 
chamber  A  served  the  dual  function  of  equilibration  vessel  and  liquid-phase 
sample  bulb.  Note  that  the  walls  of  the  liquid-phase  sample  bulb  in  Fig.  1  are 
exposed  only  to  the  liquid  phase. 

Prior  to  the  preliminary  measurements  reported  in  1976,  initial  tests  were 
conducted  to  estimate  appropriate  equilibration  periods.  In  carrying  them  out 
it  was  assumed  that  the  longest  times  would  be  required  at  low  temperatures 
because  of  the  higher  equilibrium  concentrations  and  the  (presumably) 
slower  kinetics  of  exchange  and  mixing.  After  20  h  seemed  to  be  roughly 
satisfactory  at  low  temperatures,  that  period  was  used  for  most  of  the  pre- 
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liminary  measurements.  The  results  for  oxygen  were  that  near  0°C  the 
solution  was  nearly  saturated  (A  was  only  0. 1 1 %  high),  but  as  the  temperature 
was  raised  the  solution  became  systematically  more  undersaturated  so  that  A 
was  0.76%  high  at  35CC.  Despite  this  systematic  error,  the  primary  con¬ 
clusions  from  the  preliminary  work,  i.e.,  the  functional  dependence  of  k  on 
T  and  the  relationships  between  solution  and  molecular  parameters,  remain 
valid. 


APPENDIX  B 

For  comparison,  the  notation  here  follows  Clarke  and  Glew.(a3)  At  any 
given  temperature  the  equilibrium  constant  is  related  to  the  standard  thermo, 
dynamic  function  changes  by 

*m  =  A 5#-^  (B-l) 

where  A  A  S£,  and  A  Hr  are  functions  of  temperature  only.  If  we  now 
assume  that  the  standard  thermodynamic-function  changes  are  well  behaved 
functions  of  the  reciprocal  of  the  absolute  temperature  \/T,  then  A Hj  can  be 
expanded  in  a  Taylor  series  about  the  reference  reciprocal  temperature  1/0. 
Then  defining 

fi  -  f  -  J  (B-2) 


so  that  /3  =  0  when  T  —  8,  and  also  denoting  the  derivatives  evaluated  at 
T  =  0by 


D,  = 


dl  A Hf] 

d/8‘  J,.0 


we  have  the  Taylor  expansion 


A//£  =  \Hg  +  Dd 3  +  \D2p2  +  •  ■  •  (B-3) 


We  also  have  the  relations 


and 


4Q  = 


F(A//^) 
dT  JP 


(B-4) 

(B-5) 

(B-6) 


.Jt  ,*w  . 
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so  that 


AC“r  =  [-/>,  +  d2(  i/<?)  +  •  •  ]i/r2  +  [-  d2  +  ■  •  ]i/r3  +  •  •  •  (B-7) 


Using  the  relation  between  A S%  and  A C£T,  we  have 


ASt  =  AAo 


ar 


(B-8) 


and  then  combining  Eqs.  (B-7)  and  (B-8)  with  (B-l)  and  (B-3),  we  have 


\  S°  4-  —  — —  + 

0  2  82  6  63 


R  In  Kv 


This  may  be  rewritten  as 

R  In  Kp  =  ar0  +  aJT  +  aJT2  +  a3/T3  + 


'  3Z>!  3  D2 

1 

5D2 

' 

+ 

2  +  26  + 

J*2  + 

_  _+. 

-  AH0  V  -  '  i 

+  g  e2  + 


J-3  T 


T 

( B-9) 

(B-10) 


which  corresponds  to  Eq.  (2).  From  Eqs.  (B-l),  (B-4),  (B-5),  (B-6),  and  (B-10) 
it  is  evident  that  the  thermodynamic-function  changes  ACJ,  A Hy,  AS^,  and 
ACpr  are  very  simply  related  to  the  coefficients  obtained  from  fitting  the 
experimental  data  to  Eq.  ( B-10)  for  the  temperature  dependence  of  equilibrium 
constants. 
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The  concentration  and  isotopic  fractionation  of  gases 
dissolved  in  freshwater  in  equilibrium  with  the 
atmosphere.  1.  Oxygen 

Brut  e  B.  Benson  and  Daniel  Krause,  Jr. 

I)|  -partim-nt  111  Physics.  Amherst  CdIIckc,  Amherst,  Massachusetts  01002 
Abstract 

Kqiialmm  anil  tallies  are  presented  Irinii  which  accurate  values  can  he  obtained,  in  any  nl 
several  sets  nl  units,  hir  the  concentration  of  oxygen  dissolved  in  freshwater  in  equilibrium 
with  the  atmosphere  from  O'  to  40°C  and  0.5  to  1. 1  atm.  files  are  I  last'd  on  s'ulucs  lor  the 
Henry  coefficient  ol  oxygen,  k„.  which  have  an  estimated  accuracy  of  0,02(/! .  Equations  are 
derived  which  relate  ku  to  equilibrium  concentrations  in  natural  waters.  The  equations  in¬ 
clude  corrections  lor  molecular  interactions  in  the  sailor  phase.  Uncertainty  about  the  best 
way  to  correct  lor  these  interactions  limits  the  estimated  accuracy  of  the  derived  values  to 
about  *  O  tIT'i  at  ()  ( .  and  0,04ci  at  10  (1,  hut  the  new  results  are  much  more  accurate  than 
salues  Irom  the  UNESCO  tallies.  Within  their  random  errors,  previous  measurements  arrrcc 
scry  well  with  the  now  results.  Under  equilibrium  conditions,  anti  betsveen  0°  and  fiO'X ] ,  the 
per  mil  dilference  betsveen  the  l,( ).:'■( )_,  abundance  ratio  in  the  dissolved  gas  and  the  air  is 
given  by  Is  =  -0.750  ■  (427/ 7’l,  where  T  is  in  kels  in  ami  the  standard  deviation  is  <0.02%<. 


A (.ciirate  values  lor  the  solubility  of 
oxygen  are  necessary  to  interpret  dis¬ 
solved  oxygen  concentrations  in  natural 
waters.  For  many  years,  values  based  on 
the  measurements  of  Winkler  (1889, 
1891)  and  Fox  ( 1909)  were  list'd  as  stan¬ 
dards.  After  the  determinations  of  Trues- 
dale  et  al.  (1955)  raised  doubts  about  the 
Winkler  and  Fox  standards,  new  mea¬ 
surements  were  undertaken  by  many 
workers,  blots  and  Benson  (1963//), 
Montgomery  et  al.  (1964),  Carpenter 
(1966),  and  Murray  and  Riley  (1969),  us¬ 
ing  different  methods,  obtained  agree¬ 
ment  within  the  scatter  of  their  data,  i.e. 
about  2:0. 2%,  for  the  solubility  of  oxygen 
in  pure  water  over  the  range  0°-35°C. 
Until  recently  these  have  constituted  the 
most  accurate  solubility  data  available  for 
any  gas. 

Results  from  solubility  measurements 
have  been  smoothed  in  many  different 
ways,  by  graphical  and  analytical  meth¬ 
ods  and  combinations  of  the  two  (see 
Benson  and  Krause  1976).  Valentiner 
(1927)  assumed  that  the  change  in  partial 
molal  heat  capacity  of  a  gas  during  solu¬ 
tion  was  independent  of  temperature  and 
derived  from  the  van't  Hoff  equation  the 
expression 

In  L  =  .4,,  +  (A,/T)  4  A.ln  T  (1) 
for  the  variation  of  the  Ostwald  coeffi¬ 


cient  with  temperature.  (Throughout  this 
paper  we  shall  use  T  for  kelvin  temper¬ 
atures  and  t  for  Celsius  temperatures.) 

Clarke  and  Clew  (1966)  presented  a 
systematic  treatment  of  the  generalized 
Valentiner  function 

In  k  =  A„  +  (A, IT)  +  Adn  T 

+  A, 7’  +  A +  .  .  .  (2) 

for  equilibrium  constants  in  general. 
Weiss  ( 1970)  used  the  oxygen  data  of  Car¬ 
penter  (1966)  and  Murray  and  Riley 
( 1969)  for  both  pure  water  and  seawater 
to  fit  series  of  this  form  in  temperature, 
with  additional  terms  for  the  effect  of  sa¬ 
linity,  to  Bunsen  coefficients  and  various 
concentration  values.  Postma  etal.  (1976) 
used  the  Weiss  results  to  calculate  oxy¬ 
gen  solubilities  vs.  temperature  and  sa¬ 
linity  in  the  UNESCO  International 
Oceanographic  Tables,  Vol.  2.  We  (Ben¬ 
son  and  Krause  1976)  have  proposed 

In  k  =  //„  +  (a, IT)  +  (aJT1) 

+  (a, IT')  +  .  .  .  "  (3) 

for  the  temperature  dependence  of  the 
Henry  coefficient  of  gases  dissolved  in 
pure  water,  and  we  found  that  use  of  the 
three-term  version  of  this  expression  re¬ 
veals  relationships  between  the  molecu¬ 
lar  properties  of  gases  and  the  thermo¬ 
dynamic  properties  of  their  solutions. 

Benson  et  al.  ( 1979)  have  obtained  very 


Oxyucn  in  fri'shiLutrr 


663 


precise  and  accurate  new  experimental 
results  for  tile  Henry  coefficient  of  oxy¬ 
gen  in  pure  water  from  0°  to  60  C.  I  lie 
random  error  in  these  measurements  is 
about  0.02</f,  which  is  an  order  of  mag¬ 
nitude  lower  than  in  prex  ions  w  ork 
Agreement  within  the  random  errors  in 
the  earlier  studies,  and  extensive  tests, 
indicate  that  the  new  measurements 
probably  are  free  from  systematic  errors 
and  that  the  accuracy  is  limited  by  tin- 
precision.  The  new  oxygen  data  show 
that  Fq.  3  prov  ides  a  much  better  repre¬ 
sentation  of  the  variation  of  k  with  T  than 
the  generalized  Valentiner  function  or 
any  other  function  that  has  been  consid¬ 
ered.  For  37  data  points  only  a  three-term 
power  series  in  17  is  required  to  give  a 
precision  of  fit  that  is  better  than  0.027? 
for  the  full  60°C  temperature  range.  Fur¬ 
thermore,  Henson  et  al.  have  shown  that 
Fq.  3  follows  from  a  systematic  thermo¬ 
dynamic  ana l  \  sis  that  is  cquix  alent  to  the 
Clarke  and  (ilew  i  1  C)(S(S i  treatment  of  the 
generali/ed  V  •h  ntiner  function. 

Our  purpose  in  this  paper  is  to  derive 
from  these  new  re  nits  very  accurate  val¬ 
ues  for  the  various  quantities  commonlv 
used  to  speeifv  equilibrium  concentra¬ 
tions  of  oxvgen  in  freshwater  from  ()  to 
40  (.'.  Kquafions  and  tables  are  presented 
which  constitute  significant  improve¬ 
ments  over  the  CNFSCO  tables.  In  ad¬ 
dition,  we  express  the  results  of  our  mea¬ 
surements  of  the  isotopic  fractionation  of 
oxvgen  during  solution  in  a  form  most 
useful  for  geochemical  studies. 

We  have  benefited  from  many  discus¬ 
sions  of  solubility  problems  with  1).  K. 
Carritt  and  J.  II.  Carpenter.  We  thank  I). 
K.  Carritt  and  H.  Battino  for  reading  the 
manuscript  critically.  Comments  from  C. 
If.  Mortimer  and  H.  F.  Weiss  have  been 
very  useful.  This  work  was  supported  in 
part  by  National  Science  Foundation 
grant  ()( IF  78-OS  1 20  and  in  part  by  ( )tfice 
of  Nav  al  Research  contract  N000I  4-78-C- 
0198. 

I'.tliiililninni  rimci  ntniliim  and 
tlir  limn/  corffirirnt 

Henry’s  law  states  that./,  -  k,x,.  where 
I,  is  the  lugacitv  of  the  ith  gas  in  the  v  apor 


ph  ase  and  v,  is  its  mole  traction  in  the 
liquid  phase.  Although  the  Henry  coef¬ 
ficient  k,  is  more  fundamental  because  of 
its  direct  relationship  to  the  changes  in 
the  chemical  potential  and  the  other  par¬ 
tial  molal  thermodynamic  functions,  oth¬ 
er  w  av  s  of  expressing  gas  solubilities  are 
more  useful  for  studies  of  natural  waters. 

Suppose  that  a  body  of  pure  water  is  in 
thermodynamic  equilibrium  with  an  at¬ 
mosphere  of  standard  composition  satu¬ 
rated  with  water  vapor  at  a  total  pressure 
/’.  For  oxvgen  we  may  write  Henry’s  law 


./  1)  *  I) 

i>ir  +  a, „i  +  ’>ii.i 


(4) 


where  n„  and  n,M  are  the  numbers  of 
moles  of  water  and  of  dissolved  oxvgen 
in  some  volume  \ o f  solution,  and  nKll  is 
the  number  of  moles  of  all  other  dis¬ 
solved  gases,  i.e.  essentially  nitrogen.  11 
t\,  and  rK  are  the  partial  molal  volumes 
of  oxvgen  and  of  the  residual  gas.  and  M 
and  f>  are  the  gram  molecular  mass  and 
density  ol  pure  water,  the  volume  V,.  at 
the  temperature  of  equilibrium  is  given 
by 


ir, .u,,,, 


I'  It  II  Hit 


f> 

M 


(5) 


Then,  solving  for  ii,,.i  in  the  numerator 
of  Fq.  4,  div  iding  the  result  by  V,.  and 
substituting  Fq.  5,  we  obtain  the  expres¬ 
sion 


>‘u  +  n, ,■/  *  Hit, i 

-I  if, ^  . 

(6) 

for  the  equilibrium  concentration  of  dis¬ 
solved  oxygen  in  moles  of  oxvgen  per 
volume  of  solution.  11  wo  assume  that 
the  Lewis  and  Randall  rule  applies  to 
this  system. 


k„.\t 


.In  —  XnJ  o. 


(7) 


where  is  the  lugacitv  pure  oxygen 
would  have  at  the  equilibrium  tempera¬ 
ture  if  its  pressure  were  the  same  as  tin- 
total  pressure  /’  ol  the  mixed  v  apor  phase. 
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Benson  et  al.  (1979)  showed  that 

/'„  =  PZ- „  (8) 

where  Z'„  is  the  eompressihility  laetor  for 
pure  oxygen  at  the  equilibrium  temper¬ 
ature  and  the  total  pressure  P.  (Benson  et 
al.  referred  to and  Z'„  as /*  and  Z*  (anti 
P  as  P,).  The  notation  is  changer!  here  be¬ 
cause,  following  geochemical  custom, 
the  asterisk  is  used  in  conjunction  with 
to  denote  something  quite  different. | 
Furthermore,  by  definition, 

P„  =  x,f.  (9) 

If  Eq.  7,  8,  and  9  are  combined,  and 
(1  -  fl,,P)  is  substituted  toi  Z',„  the  result 
is 

/„  =  P„(  1  -  0„P)  ( 10) 

where  is  the  negative  of  the  seconti 
pressure  coefficient  in  the  virial  expan¬ 
sion  for  the  real  gas  behavior  of  oxygen. 
The  argument  leading  to  Eq.  10  is  essen¬ 
tially  equivalent  to  assuming  that  for  a 
mixture  of  real  gases,  the  virial  cross-term 
Bu  is  equal  to  (B„  +  Bu)t 2,  where  the  sec¬ 
onti  virial  coefficient  B„  takes  account  of 
interactions  between  pairs  of  molecules 
of  type  i,  Pj,  between  pairs  of  type.;,  and 
Bh  between  mixed  pairs  of  molecules 
(e.g.  sec  Guggenheim  1950,  Sect.  5.14 
and  6. 14). 

When  Eq.  10  is  inserted  into  Eq.  6,  and 
it  is  recognized  that  the  ratio  in  the  brack¬ 
ets  in  Eq.  6  differs  from  unity  by  onl  JO 
or  30  ppm  (see  Enns  et  al.  1965  for  values 
of  the  partial  mola!  volumes),  the  result 
is 

=  "V"  =  ,,r"  (1  »„P).  (ID 

V  s  ****** 

Therefore,  at  any  P,,  and  P  the  concen¬ 
tration  of ’dissolved  oxygen  can  be  calcu¬ 
lated  from  values  for  k,„  (/,„  and  p.  from 
Benson  et  al.,  k„  is  given  by 

In  k„  -  3.71814  +  (5596.17/7  ) 

( 1.049.668'/-)  (12) 

where  T  is  in  kclvin  and  k„  m  atm.  and 
ford  -  I  -  4<rC. 

tt„  0.000975  -  (1.426  x  10  d) 

»  (6.436  x  10  »/'-).  (13) 


(E<|uation  13  is  eq.  10  in  Benson  et  al., 
which  is  bast'd  on  second  virial  coeffi¬ 
cients  given  by  Sengers  et  al.  1972.) 
Furthermore,  we  have  found  that  over 
the  range  0°-40°C  the  density  of  water 
in  g  cm  can  be  calculated  from 

In  p  =  -0.589581  +  (326.785 IT) 

-  (45,284. 1/T2)  (14) 

with  an  error  no  larger  than  20  ppm. 
Again,  T  is  in  kelvin. 

It  should  be  noted  that  Eq.  1 1  is  not 
restricted  to  the  equilibrium  of  air  with 
water.  Within  the  stated  assumptions  it  is 
valid  for  any  mixture  including  pure  oxy¬ 
gen. 

The  Bunsen  eoefheient 

A  Bunsen  coefficient  for  oxygen,  f3u, 
may  be  defined  as  the  number  of  moles 
of  dissolved  oxygen  (or,  equivalently,  the 
volume  the  dissolved  oxygen  would  oc¬ 
cupy  at  STB  if  it  were  an  ideal  gas)  per 
unit  volume  of; solution  per  unit  partial 
pressure  of  oxygen  above  the  solution. 
Various  workers  have  defined  Bunsen 
coefficients  in  other  ways,  but  this  for¬ 
mulation  is  most  useful  for  geochemical 
work.  From  this  definition  and  Eq.  If. 

=  (15i 

This  shows  that  depends  upon  not 
only  the  temperature  through  k,„  p,  and 
tfn,  but  also  the  total  pressure  P.  Usually, 
the  Bunsen  coefficient  is  referred  to  a  sit¬ 
uation  in  which  pure  oxygen  with  a  par¬ 
tial  pressure  of  1  atm  is  above  the  solu¬ 
tion,  but  Eq.  15  implies  that  the  same 
result  would  be  obtained  with  1  atm  pres¬ 
sure  of  air.  Then,  denoting  this  Bunsen 
coefficient  by  fi\„ 

d1,.  -  k  j  I  -  Poll  *  I  16) 

where  /’„.,  is  the  saturated  v  apor  pressure 
of  water  in  atmospheres  at  the  tempera¬ 
ture  of  equilibrium. 

Because  the  present  discussion  is  lim¬ 
ited  to  temperatures  below  40'C.  the 
quantitv  can  be  neglected.  (Even 

at  40  C,  it  is  0.00004.!  Eq.  16  then  re¬ 
duces  to 
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/3'" =  k\i(l  ~ w,,)'  (17> 

Although  the  inn. lining  factor  ( 1  -  9,,) 
differs  onl v  slightly  from  unity  tit  varies 
from  0.99902  at  0°C  to  0.99949  at  40°C), 
for  accurate  values  of /E)1,,  it  must  he  in¬ 
cluded  to  correct  for  molecular  interac¬ 
tions  in  the  vapor  phase.  Substituting 
M  =  18.0153  g- mol  '  in  Eq.  17  gives 

0\.  =  5.5508  x  10-/'-  (1  -  ttj  (18) 
*(> 

where /S', ,  is  inpmoldm  '*  of  solution  per 
atm  when  p ,  the  density  of  pure  water,  is 
in  g  em  3  and  k„  is  in  atm.  For  /31,,  in  cm'* 
of  ideal  gas  at  STP  dm  3 -atm  we  re¬ 
place  the  numerical  coefficient  in  Eq.  18 
with  1.24410  x  10'*.  (The  volume  of  1 
mole  ol  ideal  gas  at  STP  has  heen  taken 
to  be  22,4  14  dll'1.) 

To  facilitate  the  determination  of  val¬ 
ues  for  /3\,,  we  have  fitted  the  logarithms 
of  values  calculated  from  Eq.  18  (with  the 
aid  of  Eq.  12-14)  to  a  three-term  power 
series  in  the  reciprocal  of  the  absolute 
temperature.  The  result  is 

lu  /3'„  =  13.5246  -  (5.26895  x  KP/T) 

+  (1.00417  x  lO'VT-)  (19) 

where  is  in  jumoldm  ■’■atm  For 
values  in  cm*  of  ideal  gas  at  STP- 
dm  ’atm  ’,  we  substitute  9.7265  for  the 
first  constant.  Over  the  range  0°— 40°C, 
values  determined  from  Eq.  19  differ 
from  those  obtained  from  Eq.  18  In  less 
th  an  tlu*  uncertainty  in  k„.  The  excellent 
fit  of  In  fj\,  to  the  three-term  series  in 
FT  is  not  surprising,  because  p  and  (1  - 
9,3  varv  only  slightly  with  temperature. 

lu/uilihrit,  in  concentrations  from  a 
unit  standard  atmosphere 

For  oxygen,  the  unit  standard  atmo¬ 
spheric  concentration  by  volume,  C  *,„  or, 
more  briefly,  the  unit  standard  concentra¬ 
tion.  is  defined  here  as  the  concentration 
of  dissolved  oxygen  per  unit  volume  of 
solution  (measured  at  equilibrium  tem¬ 
perature)  when  it  is  in  equilibrium  with 
an  atmosphere  of  standard  composition 
and  saturated  with  water  vapor,  at  a  total 
pressure  (including  that  of  the  water  sa¬ 


por)  of  1  atm.  In  the  past  solubility,  equi¬ 
librium  solubility,  and  sometimes  air  sol¬ 
ubility  have  been  used  for  this  concept. 
The  first,  however,  also  is  commonly 
used  to  refer  in  a  general  collective  sense 
to  the  solution  process  and  ways  of  de¬ 
scribing  it.  Neither  solubility  nor  equilib¬ 
rium  solubility  suggests  that  the  system 
is  composed  of  the  solution  and  the  at¬ 
mosphere,  while  air  solubility  literally 
implies  a  concentration  of  air.  Although 
no  procedure,  short  of  creating  a  new 
word,  is  ideal  for  both  clarity  and  brevity, 
choice  of  unit  standard  concentration  of¬ 
fers  the  advantages  of  indicating  that  it  is 
a  concentration  under  standard  condi¬ 
tions,  i.e.  equilibrium  with  a  standard 
atmosphere  at  unit  pressure.  The  unit 
standard  concentration  by  mass,  Ct,,,  is 
defined  in  an  analogous  way,  except  that 
the  concentration  is  referred  to  unit  mass 
of  solution.  If  both  C*,,  and  Ct„  are  re¬ 
ferred  to  as  unit  standard  concentrations, 
the  asterisk  and  dagger,  and  the  associ¬ 
ated  units,  will  distinguish  between 
them. 

In  order  to  write  down  an  expression 
for  from  Eq.  1 1,  we  need  an  expres¬ 
sion  for  P,,.  For  temperatures  below  40°C 
the  saturated  vapor  pressure  of  pure 
water  is  <0.08  atm.  Consequently,  the 
common  practice  of  taking  the  difference 
between  the  total  pressure  P  =  1  atm  and 
the  value  for  P,,.,.  read  from  the  steam  ta¬ 
bles  (in  atm)  probably  yields  a  good  ap¬ 
proximation  for  the  partial  pressure  of  the 
remaining  gases.  Then,  again  using  the 
definition  that  the  partial  pressure  of  a 
gas  in  a  mixture  is  equal  to  its  mole  frac¬ 
tion  times  the  pressure  of  the  mixture,  we 
have 

P„  =  0.20946(  1  -  P„,.),  (20) 

and 

<•*„  -  0.20946(1  -  P„.,.)  ( I  -  9„). 

K{tM 

(21) 

or,  in  terms  of  from  Eq.  17, 

C*„  =  0.20946(1  -  P, 3',,.  (22) 

The  xalue  0.20946  for  the  mole  fraction 
of  oxygen  in  dry  air  is  from  Cllueckauf 
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(1951).  For  convenience  in  call  ulating 
we  have*  expressed  the  temperature 
dependence  of  bv 

In  =  11.8571  -  (.1,840. 70/7') 

(210,901/7'-)  (23) 

where  T  is  in  kelvin  and  is  in  atm. 
Between  (F  and  4(FC  Eq.  23  gives  value's 
for  ( 1  )  that  differ  I'roin  the  standard 

values  l>y  truly  I  1  ppm  rms.  Using  fc„,  ft,,, 
p,  and/’,,,  from  Kq.  12,  13,  14,  and  23,  we 
have  ealeulated  values  for  from  Eq. 
21  at  integral  temperatures.  These  are 
given  in  Table  1  in  three  sets  of  units: 
pg-atoms  dm  ",  cm"  (ideal  gas  at  STP)- 
dm  ",  and  mg  dm  ".  |To  eonforin  with 
the  International  Association  of  Physical 
Oceanographers  ( I lclland-1  lansen  et  al. 
1948)  and  the  Joint  Panel  on  Oceano¬ 
graphic  Tables  and  Standards  (Postma  et 
al.  1970)  we  have  expressed  the  first  set 
in  pg-atoms  rather  than  pmol.  For  the  last 
set  1  mole  of  oxvgcn  has  been  taken  to 
be  31.9988  g.| 

To  obtain  it  is  necessary  only  to 

divide  (■’*,,  In'  the  density  of  the  solution. 
To  an  excellent  approximation,  however, 
the  latter  is  equal  to  p,  the  density  of  pure 
water.  Therefore,  from  F,q.  21, 

C* 

ft  =  11 

P 

=  0.20940(1  -  1  (1  -  »„). 

k  i » • » * 

(24) 

Values  for  have  been  calculated  at 
integral  temperatures  from  Eq.  24.  They 
are  shown  in  the  last  three  columns  of 
Table  I  in  pg-atoms  kg  ',  cm"  (ideal  gas 
at  STP) -kg  '.and  mgkg  '. 

For  use  in  either  fieldwork  or  the  lab¬ 
oratory  Table  I  provides  a  simple  means 
for  finding  very  accurate  values  for 
and  Linear  interpolation  between 

degrees  introduces  an  error  not  >0.017%. 
For  computer  application  we  have  fitted 
both  In  and  In  calculated  as 

above,  to  power  series  in  the  reciprocal 
of  the  kclvin  temperature.  Although,  as 
we  have  seen.  In  /71,,  fits  a  three-term 
function  extraordinarily  well,  more  terms 
are  required  for  and  C.'t,,  because  (1 


/’„,)  has  a  different  functional  form 
from  The  standard  deviations  for 

three-,  lour-,  and  five-term  series  are 
0.16,  0.026,  and  0.0047, .  The  two  five- 
term  series  are  the  briefest  within  the 
precision  of  the  Benson  et  al.  (1979)  val¬ 
ues  fork,,.  For  in  pg-atoms  •  dm  ", 

In  C*„  =  -  135.20890 

+  (1.575701  x  lOVT) 

(6.642308  x  107/T’) 

+  ( 1.2438(H)  x  10"V7'") 
-(8.621949  x  10"/7’4).  (25) 

To  obtain  in  cm"- dm  1  or  mg- dm  ", 
we  replace  the  first  constant  in  Eq.  25 
with  -  139.7001 1  or  -139.34410.  Similar¬ 
ly,  forCt,,  in  pg-atoms- kg 

In  C1\,  =  -134.60690 

+  (1.572288  x  10s//') 

-  (6.637149  x  1()7P) 

+  (1.243678  x  10"'/T") 

-  (8.621061  x  1()"/T4),  (26) 

and  when  C.'t,,  is  desired  in  cm"- kg  1  or 
mg  kg  ',  we  substitute  -139.09811  or 
-  138.74210  for  the  first  constant  in  Eq.  26. 

Equilibrium  concentrations  vs. 
atmospheric  pressure 

It  is  useful  to  have  a  simple  way  of  cal¬ 
culating  the  equilibrium  concentration  of 
dissolved  oxygen  in  a  situation  where  the 
atmospheric  pressure  is  different  from 
unity.  The  more  general  form  of  Eq.  21. 

(:'•„  =  0.20946(P  -  />„,.)  /  (1  -  ftt>P). 

K‘0.M 

(27) 

is  the  appropriate  equation  to  use  (if  we 
assume  that  the  composition  of  the  at¬ 
mosphere  is  uniform),  but  to  use  it  in 
the  field  would  be  cumbersome.  How¬ 
ever, 

rr  =  r*  pH1  -  P,n-/P)(\  -  0„P) 

'  "  '"l  (I  -  JVrKl  -  #„')  I' 

(28) 

Values  for  the  quantity  in  brackets  for 
oxygen  are  given  in  Table  2  for  pressures 
from  1.1  to  0.5  atm  (the  latter  corresponds 
to  an  altitude  of  about  5  km),  and  for  tem¬ 
peratures  from  0°  to  4()°C.  Accurate  val- 


Oxygen  in  freshwater 


667 


Table  I  Willies  lor  (  *,,  and  (  +,,  in  various  sets  ol  units  at  integral  temperatures.  Linear  interpolation 
between  decrees  will  introdnee  an  error  •  O.OIT's.  To  eonvert  from  pit-atoms  to  prnol  div  ide  by  2. 


(  *.. 

lemp 

mi-.dmus 

mil 

p  il-.iloms 

mg 

C 

tin)1 

•w 

377 

kg 

!g 

kg 

0.0 

913.85 

10.2415 

1  1.021  1 

913.99 

10.2432 

14.0235 

1.0 

S8S.53 

9.0577 

11.2100 

SSS.02 

9.9588 

14.2174 

2.0 

804.30 

9. 8809 

13.829.3 

SOI.  41 

9.0875 

13.8302 

1.0 

84 1 ,2S 

9.4282 

1.1.40(H) 

SI  1.31 

9.1280 

13.4000 

TO 

S  10.2.1 

9.181  1 

13.1072 

S  19.20 

9.1814 

13.1077 

5.0 

70S.  1.5 

8.9448 

12.7090 

70S.  IS 

8.9452 

12.7705 

0.0 

777  OS 

8.7188 

12.4173 

778.03 

S.719-1 

12.4481 

7.0 

75S.0S 

8.5025 

12.1385 

75S.75 

S. 50.34 

12.1397 

SO 

740  10 

8.29.5-1 

1  1  .8427 

740.31 

S.2900 

1 1 .8445 

0.0 

722. 4S 

8.0i)89 

1 1.5593 

722.04 

S.09S0 

1 1.5019 

10.0 

705.50 

74)085 

1 1.2870 

705.7 1 

7.9089 

1 1.2910 

I  1.0 

(3S9  2  1 

7.7240 

1 1.0270 

089.48 

7.7270 

1 1.0.313 

12.0 

074  57 

7.5487 

I0.770S 

073.91 

7,5525 

10.7822 

no 

058.50 

7.4805 

10.5300 

0583X1 

7.3850 

10.543 1 

1  1.0 

044  1  1 

7  2187 

10  30.57 

044.01 

7.2242 

10.3134 

15.0 

0  10.25 

7.0844 

10.0S.37 

0.30.82 

7.0090 

10.0927 

io.o 

010,01 

8.9147 

9.S702 

017.55 

0.9210 

9.8800 

17.0 

004. IMi 

8.7897 

9.00-17 

004.80 

0.7780 

9.0704 

IS.O 

50 1  00 

8  841  1 

9.1007 

592.51 

0.0403 

9.4799 

10.0 

.570.70 

0.4974 

9.2759 

580.08 

0.5078 

9.2907 

20.0 

508.27 

8  1888 

9.0920 

509.28 

0.38(H) 

9.108,3 

21.0 

557  17 

8  2444 

S.9I45 

558.29 

0.2508 

8.9.324 

22.0 

540  to 

8  1242 

8.7431 

547.08 

0.1379 

8.7020 

25.0 

540. 12 

8.0085 

8.5770 

537.44 

0.0231 

8.5987 

24.0 

520. 1 1 

8.8982 

S.  1 170 

527.54 

5.9122 

8.4404 

25.0 

510.44 

5.7878 

S.202S 

517.97 

5.8049 

8.2873 

20.0 

507.  OS 

5.0828 

S.  1 130 

508,71 

5.7012 

8.1392 

27.0 

4984)1 

5.5812 

7.9079 

499.7.5 

5.0007 

7.9958 

2S.0 

ISO  22 

5  1827 

7.S272 

491.07 

5.50.34 

7.8509 

20.0 

tso.oo 

5  4871 

7  090S 

482.0.5 

5.4091 

7.7222 

.100 

172  42 

5  2944 

7  5.584 

474.49 

5.3 1 70 

7.5915 

•41.0 

104.  is 

5.2044 

7  I29S 

400.50 

.5.2287 

7.4(347 

12.0 

150  50 

5i  187 

7.30  IS 

458.85 

5.1424 

7.3414 

1.4.0 

its  00 

5.0415 

7. 1  S3 1 

451.30 

5.0584 

7.2215 

110 

1 1 1  50 

4  9188 

7. (Hi  17 

1 14  (Hi 

4.9707 

7.1048 

15.0 

1.11.45 

1.8877 

0.949.3 

1.30.90 

4.8970 

0.9912 

10.0 

127.11 

1.7889 

0.S30S 

430.04 

4  8194 

0.8803 

17  0 

120  15 

4  7120 

0.7209 

42.3.28 

4.7437 

0.7722 

ISO 

11171 

4.8488 

0.0190 

1 10.08 

4.0097 

0.0000 

40.0 

107  I  S 

4  5844 

0.51  17 

110.22 

4.5974 

0.50.34 

to  0 

1(H)  70 

1  4914 

0.4  1  19 

40.3.91 

4.5200 

0.4023 

ties  for  tail  he  obtained  by  simply 

multiplying  from  eitlier  Table  I  or 
Eq.  25,  hy  tile  product  of  /’  and  the  ap¬ 
propriate  value  for  the  bracket  from  Ta¬ 
ble  2. 

Linear  interpolation  in  F  and  /  will  in¬ 
troduce  an  error  •  0.02*7  in  the  upper  and 
left  sections  of  Table  2.  Interpolations 
from  numbers  in  parentheses  will  lead  to 
errors  - 0.05 Of.  With  the  numbers  in 
brackets,  the  interpolation  errors  rapidly 


become  larger,  but  this  region  for  P  and 
t  is  not  usually  of  interest  in  geochemical 
studies.  It  is  obvious  that  Eq.  28  and  the 
procedure  just  described  are  equally  ap¬ 
plicable  to 

Isotopic  fractionation  of  atmospheric 
oxygen  (luring  solution 

From  direct  mass  speetrometric  com¬ 
parison  of  the  ’’O.iO;  abundance  ratio 
in  the  dissolved  gas  with  the  eorrespond- 
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Table  2.  Values  lor  quantity  bracketed  in  Kq.  28  \s.  V  and  t .  Linear  interpulatum  in  /'  and  t  will 
introduce  an  error  ■  0.02'/!  in  upper  and  left  sections  of  table  Interpolations  using  numbers  in  paren¬ 
theses  will  lead  to  errors  - "0.059S .  With  the  numbers  in  brackets,  inter|><>lation  errors  In-c  ome  larger 


Either  temperature  or  pressure  may  he  interpolated  first,  as 
by  the  two  arrays  below'  tin*  body  of  the  table. 

illustrated  lor 

1  IDO  L  and  V 

0.67  atm 

P.  atm 

•t: 

11 

1.0 

0« 

us 

0.7 

Of) 

0  5 

0.0 

1.0005 

l.(MHK) 

0.9994 

0.9987 

0.9977 

0.9963 

0.9944 

5.0 

1.0007 

1.0000 

0.9991 

0.9980 

0.9966 

0.9946 

0.9918 

10.0 

1.0010 

l.(XXX) 

0.9957 

0.997 1 

0.9950 

0.9922 

0.9882 

15.0 

1.0015 

l.(XXX) 

0.9982 

0.9959 

0.9929 

(0.9889) 

(0.9833) 

20.0 

1 .002 1 

l.(KKH) 

0.9974 

0.9942 

(0.9901) 

(0.9845) 

10.9767) 

25.0 

1.1X129 

l.(XMK) 

0.9965 

0.9921 

(0.9864) 

(0.9787) 

|0.9680| 

,10.0 

1 .0019 

1.0000 

0.9952 

(0.9892) 

(0.9814) 

10.9711) 

10.9566) 

45.0 

1 .0053 

1.0000 

(0.9935) 

(0.9854) 

(0.97.50) 

10.96101 

(0.94151 

40.0 

1.0071 

l.(XXX) 

(0.9913) 

(0.9805) 

10.9665) 

10.9479) 

10.9217) 

Temp  first 

Prt*ss  first 

0.7 

0.ft7 

0.6 

0  7 

0.67 

0.6 

0 

3 

0.9977 

0.99704 

0.9965 

0.9963 

0.99528 

0 

3 

0.9977 

0.99728 

0.9965 

0.9963 

5 

0.9966 

0.9946 

5 

0.9966 

0.99600 

0.9946 

ini'  ratio  in  the  gas  above  the  solution, 
Benson  et  al.  ( 1979)  confirmed  the  earlier 
observations  by  Kbits  and  Benson  (1963/>) 
on  oxygen  in  pure  water,  and  by  Kroop- 
niek  and  Craig  (1972)  on  oxy  gen  in  sea¬ 
water,  that  at  equilibrium  the  heavier 
molecule  is  slightly  more  soluble.  In  ad¬ 
dition,  however,  the  higher  precision  and 
the  wider  temperature  range  of  the  new 
measurements  provide  quantitative  in¬ 
formation  about  the  variation  of  the  oxy¬ 
gen  isotope  fractionation  with  tempera¬ 
ture.  These  results  are  directly  applicable 
to  atmospheric  oxygen,  because  the  pres¬ 
ence  of  the  other  gases  should  not  influ- 
epee  the  fractionation  significantly,  but 
for  geochemical  purposes  the  most  useful 
way  to  express  the  fractionation  as  a  func¬ 
tion  of  temperature  is 

fi  =  -0.7.30  +  (427/7’)  (29) 

where  T  is  in  kelvin,  and  S  is  the  per  mil 
difference  between  the  34()a:S202  abun¬ 
dance  ratio  in  the  dissolved  gas  and  the 
ratio  in  atmospheric  oxygen 
when  freshwater  is  in  equilibrium  with 
air. 

We  emphasize  that  Kq.  29  will  give  the 
isotopic  difference  between  dissolved 
and  atmospheric  oxygen  under  ei/uilib- 


rium  conditions  only.  In  addition  to  bi¬ 
ological  and  chemical  activity,  physical 
kinetic  effects  can  influence  the  isotopic 
ratio  in  the  dissolved  gas.  For  example, 
during  our  preliminary  measurements 
(Benson  and  Krause  1976),  we  deter¬ 
mined  the  isotopic  ratio  in  the  dissolved 
gas  as  a  function  of  the  interval  after  the 
degassed  water  was  brought  into  contact 
with  the  gas.  Initially,  as  a  result  of  the 
faster  kinetics  of  the  lighter  molecule,  the 
:u( b./'-'t b.  ratio  in  the  dissolved  gas  was 
found  to  be  lower  than  in  the  gas  above 
the  solution,  but  as  time  progressed  the 
dissolved  ratio  increased,  then  surpassed 
that  of  the  gas  above,  and  finally  stabi¬ 
lized  at  the  higher  equilibrium  value. 

Discussion 

The  empirical  method  described  by 
Benson  et  al.  (1979)  determines  k  direct¬ 
ly.  It  yields  fully  corrected  sallies  for  k 
without  the  necessity  of  making  theoret¬ 
ical  assumptions  about  how  to  take  care 
of  molecular  interactions  in  the  vapor 
phase.  Similarly,  the  methods  used  by 
Montgomery  et  al.  (1964),  Carpenter 
(1966),  and  Murray  and  Riley  (1969)  de¬ 
termine  C*  directly.  To  calculate  C*  from 
k.  or  vice  versa,  however,  it  is  necessary 
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to  include  virial  corrections  if  the  data  are 
sufficiently  accurate.  In  the  absence  of 
experimental  information  about  cross-in¬ 
teractions  in  gas  mixtures ,  we  have  in¬ 
voked  the  Lewis  and  Randall  rule  and 
assumed  the  applicability  of  the  defini¬ 
tion  of  partial  pressure  which  is  implicit 
in  the  expressions  P0  =  xt)P  and  Pt,  = 
0.20946(P  -  P, Other  theoretical  pos¬ 
sibilities  exist.  For  example,  in  our  ex¬ 
perimental  measurements  of  k  we  initial¬ 
ly  tried  to  correct  for  interactions  between 
gas  and  water  molecules  in  the  vapor 
phase  by  calculating  the  virial  cross-term 
Btj,  using  methods  discussed  in  Hirseh- 
felder  et  al.  (1967).  The  residts  were  not 
adequate  for  the  high  precision  of  the 
measurements,  and  we  finally  adopted 
the  empirical  approach. 

In  Fig.  1  we  illustrate  the  significance 
of  these  corrections  and  compare  the  ex¬ 
perimental  values  for  C*„  from  Mont¬ 
gomery  et  al..  Carpenter,  and  Murray  and 
Riley  with  those  calculated  from  Eq.  21 
using  our  values  of  ktt.  (The  results  from 
Klots  and  Benson  (1963o)  are  not  includ¬ 
ed  here  because  vve  want  to  focus  on  the 
directly  measured  C*„  values.  Benson  et 
al.  ( 1979)  showed  that  the  results  of  Klots 
and  Benson,  Carpenter,  and  Montgomery 
et  al.  are  quite  similar.)  The  reference 
line  in  Fig.  1  is  based  on  our  smooth  val¬ 
ues  for  k„,  and  the  tight  fit  of  the  solid 
circles  to  it  indicates  the  high  precision 
of  the  measurements.  Before  plotting  the 
percentage  deviations  of  the  experimen¬ 
tal  values  forC.'*,,  from  the  three  sources, 
we  expressed  their  values  in  /ug- 
atoms  dm  '1  as  follows.  Montgomery  et 
al.  expressed  their  results  as  the  number 
of  milligrams  of  dissolved  oxygen  per  li¬ 
ter  of  solution,  with  the  volume  deter¬ 
mined  at  20"C  rather  than  at  the  temper¬ 
ature  of  equilibrium.  We  have  used 

C ==  C'"'  x  2  x  ia'  x  p  , 
31.9988  p.  ’ 
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Fig.  1.  Percentage  deviations  from  our  present 
smoothed  values  for  C*{t,  in  /ugsitoms  dm  *,  of  the 
measured  values  from  Montgomery  et  al..  Carpen¬ 
ter,  and  Murray  and  Riley,  and  of  the  UNESCO 
smoothed  values.  Dashed  lines  indicate  possible 
range  of  uncertainty  of  our  present  smoothed  val¬ 
ues. 


factor  89.28  fig-atoms  ml1.  Murray  and 
Riley  stated  that  they  used  a  modification 
of  the  method  of  Carpenter  (1965),  in 
which  the  factor  5,598  ml  of  02  per  equiv¬ 
alent  is  given.  Accordingly,  we  have  mul¬ 
tiplied  the  Murray  and  Riley  experimen¬ 
tal  sallies  (by  the  Winkler  method  only) 
by  89.318  fig-atoms -ml'1. 

The  dashed  lines  in  Fig.  1  indicate  the 
range  of  possible  virial  corrections  for 
oxy  gen.  If  no  correction  were  made,  i.e. 
if  the  factor  (1  -  #„)  in  Eq.  21  were  taken 
to  be  unity,  the  reference  line  and  the 
solid  circles  would  move  upward  (rela- 
tive  to  the  other  three  sets  of  results)  to 
the  position  of  the  upper  dashed  line. 
The  shift  would  vary  from  about  0.1 %  at 


Table  3.  Comparison  of  deviations  oi  experi¬ 
mental  values  for  C*„  from  our  present  results. 


where  p  and  p..„  are  the  densities  of  water 
at  the  equilibrium  temperature  and  20°C, 
to  obtain  the  corresponding  values  tor 
C*,,.  Only  values  from  Carpenter’s  table 
1  have  been  included  here,  and  they 
have  been  multiplied  In  bis  conversion 
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0°C  to  0.05%  at  40°C.  On  the  other  hand, 
if  we  assumed  that  the  correction  for  the 
air-water  system  were  the  same  as  the 
empirical  correction  found  hy  Benson  et 
al.  (1979)  for  the  pure  oxygen-water  sys¬ 
tem,  then  (1  -  0„)  would  he  replaced  by 
[1  -  (A, i  -  0„)],  where  A()  is  given  by  Ben¬ 
son  et  al.  as  a  function  of  temperature. 
With  this  assumption  the  reference  line 
would  become  the  lower  dashed  line, 
which  is  about  0.085%  lower  at  0°C  and 
0.03%  lower  at  40°C. 

In  principle,  comparison  of  the  results 
from  the  two  kinds  of  measurements 
could  test  the  possibilities  discussed 
above,  but  it  is  clear  from  Fig.  1  that  the 
scatter  in  the  three  sets  of  previous  data 
is  too  large.  A  precision  of  about  0.02% 
for  the  direct  C*„  measurements  would 
be  required  for  such  a  test.  Nevertheless, 
within  their  random  errors  the  older  data 
are  in  general  agreement  with  our  re¬ 
sults.  The  Montgomery  et  al.  and  Car¬ 
penter  values  are  distributed  approxi¬ 
mately  equally  above  and  below  the 
reference  line,  while  the  Murray  and  Bi- 
ley  points  tend  to  be  slightly  lower. 
These  comparisons  are  expressed  in 
another  way  in  Table  3,  which  shows  that 
the  Montgomery  et  al.  and  Carpenter  de¬ 
viations  each  average  essentially  to  zero 
with  an  rms  deviation  of  about  0.2%,  and 
the  Murray  and  Hiley  deviations  average 
0.15%  low  with  an  rms  deviation  of 
0.25%.  Table  3  differs  somewhat  from 
Table  3  in  Benson  et  al.,  primarily  be¬ 
cause  there  the  |1  -  (A,,  -  0J|  correction 
was  applied  to  the  three  sets  of  experi¬ 
mental  (.'*,>  values,  and  a  different  con¬ 
version  factor  was  used  for  the  Murray 
and  Riley  values. 

Although  we  cannot  determine  the 
best  correction  to  choose  for  calculating 
C*,,  from  values  for  k„,  the  two  dashed 
lines  in  Fig.  1  provide  probable  limits  of 
uncertainty  and  permit  an  estimation  of 
the  accuracy  of  the  values  in  Table  1.  Be¬ 
cause  the  second  virial  coefficient  for  ni¬ 
trogen  is  smaller  than  that  for  oxygen,  the 
[1  -  (A,,  -  0,,)|  correction  probably  is  too 
large.  On  the  other  side,  it  is  improbable 
that  no  virial  correction  would  be  re¬ 
quired.  Including  the  possible  0.02%  sys¬ 
tematic  error  in  k„  given  by  Benson  et  al.. 


we  estimate  roughly  that  the  values  of 
C  and  Ct(,  in  Table  1  are  good  to  better 
than  ±0.07%’  at  lower  temperatures  and 
perhaps  ±0.04%  near  40°C.  These  uncer¬ 
tainties  are  negligible  for  studies  of  nat¬ 
ural  waters. 

Weiss  (1970)  fitted  the  Carpenter 
(1966)  and  Murray  and  Riley  (1969)  data 
for  both  freshwater  and  seawater  to  four- 
term  Valentiner-type  functions  for  the 
variation  of  and  with  tempera¬ 
ture,  with  an  additional  three  terms  for 
the  dependence  on  salinity.  Mis  values 
lie  below  our  results  by  only  0,05-0.10% 
at  low  temperatures.  Although  they  de¬ 
viate  more  significantly  above  20°C,  the 
values  of  Weiss  still  agree  with  ours  with¬ 
in  the  uncertainties  of  the  previous  data. 

In  1976,  Postma  et  al.  constructed  the 
UNESCO  International  Dceanonrapliie 
Tables  from  the  equations  of  Weiss.  Un¬ 
fortunately,  they  apparently  assumed  that 
the  values  of  Weiss  were  given  as  the  vol¬ 
ume  the  gas  would  occupy  at  STP  if  it 
were  ideal.  (The  factor  89.23  used  by 
Postma  et  al.  to  convert  from  volume  of 
gas  at  STP  to  /ag-atoms  is  equivalent  to 
22,414  cm'’- mol  '.)  In  fact,  Weiss  (pers. 
comm.)  intended  his  values  to  be  ex¬ 
pressed  in  terms  of  the*  volume  which 
would  be  occupied  by  a  real  gas  at  STP. 
(The  corresponding  conversion  factor 
would  be  89.318,  which  is  equivalent  to 
22,392  cnr’  ino!  >.)  The  UNESCO  values 
are  systematically  low  by  from  0.15  to 
0.35%.  This  is  shown  in  Fig.  1,  which  in¬ 
cludes  a  plot  of  the  deviations  of  the 
UNESCO  values  from  ours. 

It  should  be  noted  that  it  is  equally 
proper  to  express  the  amount  of  dissolved 
gas  in  terms  of  its  STP  volume  as  a  real 
gas  or  as  a  hypothetical  ideal  gas,  provid¬ 
ed  that  it  is  made  clear  which  is  being 
nstal.  We  prefer  the  ideal  gas  notation, 
because  the  values  specified  do  not  de¬ 
pend  on  the  state  of  our  knowledge  of  the 
virial  coefficients  of  the  gas.  The  best 
procedure,  however,  is  to  eliminate  the 
source  of  confusion  by  expressing  the 
amount  of  dissolved  gas  directly  in  /u.g- 
atoms  (or  jxmol)  or  mg.  Other  reasons  for 
this  choice  are  discussed  by  Helland- 
Ilansen  et  al.  (1948)  and  Carpenter 
(1966),  Wo  have  included  the  volume 
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specification  in  our  equations  and  tables 
only  because  of  its  historical  use. 

We  pointed  out  earlier  that  the  Henry 
coefficient  k  is  the  quantity  most  funda¬ 
mental  to  describing  the  solution  pro¬ 
cess.  On  the  other  hand,  for  many  geo¬ 
chemical  and  engineering  purposes  C* 
(or  Ct)  is  the  directly  applicable  quantity. 
To  obtain  the  Bunsen  coefficient  from 
either  k  orC*  requires  assumptions  about 
the  behavior  of  the  vapor  phase.  Further¬ 
more,  because  of  molecular  interactions 
in  the  vapor  phase,  the  Bunsen  coeffi¬ 
cient  is  not  a  constant  at  a  given  temper¬ 
ature,  but  is  a  function  of  the  total  pres¬ 
sure  (Eq.  15).  This  means  that  for  precise 
work  the  Bunsen  coefficient  loses  its  pri¬ 
mary  usefulness,  because  we  cannot  sim¬ 
ply  multiply  it  by  the  partial  pressure  of 
the  gas  above  a  solution  to  obtain  the 
equilibrium  concentration  of  the  gas  in 
the  solution.  In  general,  it  is  necessary  to 
use  Eq.  11  to  calculate  C  or,  when  the 
gas  mixture  is  air,  C*,  Ct,  or  Cp  can  be 
used  as  discussed  above.  For  these  rea¬ 
sons,  we  suggest  that  in  describing  the 
solubility  of  a  gas,  emphasis  should  be 
placed  on  k,  C*,  or  Ct  rather  than  on  the 
Bunsen  coefficient.  Avoiding  the  latter 
would  have  the  additional  advantages  of 
eliminating  the  confusion  associated 
with  the  many  different  ways  it  has  been 
defined  and  of  settling  the  conflict  over 
whether  to  use  the  symbol  a  for  the  Bun¬ 
sen  coefficient  or  for  the  isotopic  fraction¬ 
ation  factor.  We  agree  with  the  sugges¬ 
tion  by  Weiss  (1970)  that  a  be  used  for 
the  latter. 

References 

Benson,  B.  B.,  and  D.  Khausk.  Jh.  1976.  Empir¬ 
ical  laws  for  dilute  aqueous  solutions  of  non¬ 
polar  gases.  J.  Cliein.  Phvs.  64:  689-709. 

- ,  - .  and  M.  A.  Peterson.  1979.  The 

solubility  and  isotopic  fractionation  of  gases  in 
dilute  aqueous  solution.  1.  Oxvgen.  J.  Solution 
Chem.  8:  655-690. 

Carpenter,  J.  H.  1965.  The  Chesapeake  Bay  In¬ 
stitute  technique  for  the  Winkler  dissolved  oxy¬ 
gen  method.  I.imuol.  Oceanogr.  10:  141-143. 
- — .  1966.  New  measurements  of  oxygen  solu¬ 
bility  in  pure  and  natural  water.  I.imnol, 
Oceanogr.  11:264-277. 

Ci.ahkk,  K.  C„  AND  O.  N.  Ci.ew.  1966.  Evaluation 
of  thermodynamic  functions  from  equilibrium 
constants.  Trans.  Faraday  Soc.  62:  539-547. 
Fans.  T.,  P.  F.  Sciioi  xnper,  and  K.  I).  Brad- 


street.  1965.  Effect  ol  hydrostatic  pressure 
on  gases  dissolved  in  water.  J.  Phvs.  Chem.  69: 
389-391. 

Fox,  C.  J.  1909.  On  the  coefficients  of  absorption 
of  nitrogen  and  oxygen  in  distilled  water  and 
seawater,  and  ol  atmospheric  carbonic  acid  in 
seawater.  Trans.  Faraday  Soc.  5:  68-87. 

Cl.UECK.AUF,  E.  1951.  The  composition  of  atmo¬ 
spheric  air.  p.  3-10.  In  Compendium  of  mete¬ 
orology.  Am.  Meteorol.  Soc. 

Cuc.CKNHKlM,  E.  A.  1950.  Thermodynamics,  an 
advanced  treatment  for  chemists  and  physi¬ 
cists,  2nd  ed.  North-llolland. 

Hei.i.ano-Hanskn,  B.,  J.  P.  Jacobsen,  and  T.  G. 
Thompson.  1948.  Chemical  methods  and 
units.  Pub!.  Sei.  Assoc.  Oceanogr.  Phvs.  9. 

28  p. 

1 1 tRSCH FEEDER,  J  ().,  C.  F.  CURTISS,  AND  K.  B. 
Bird.  1967.  Molecular  theory  of  gases  and  liq¬ 
uids.  Wiley. 

Ki.ots,  C.  E.,  and  B.  B.  Benson.  196.3«.  Solubil¬ 
ities  of  nitrogen,  oxygen,  and  argon  in  distilled 
water.  J.  Mar.  Kes.  21:  48-57. 

- ,  AND  - .  196.3b.  Isotope  effect  in  the 

solution  of  oxygen  and  nitrogen  in  distilled 
water.  J.  Chem.  Phvs.  38:  890-892. 

Kroopnick,  P.,  and  H.  Craic.  1972.  Atmospheric 
oxygen:  Isotopic  composition  and  solubility 
fractionation.  Science  175:  54-55. 

Montgomery,  H.  a.,  N.  S.  Thom,  and  a.  Cock- 
burn.  1961.  Determination  of  dissolved  oxy¬ 
gen  by  the  Winkler  method,  and  the  solubility 
of  oxygen  in  pure  water  and  in  seawater.  J. 
Appl.  Chem.  14:  280-296. 

Murray,  C.  N.,  and  J.  1*.  Hii.ky.  1969.  The  solu¬ 
bility  of  gases  in  distilled  water  and  seawater- — 
2.  Oxygen.  Deep-Sea  Res.  16:  311—320. 

Postma,  H.,  A.  Svansson,  H.  I.acombe,  and  K. 
Grasshoff.  1976.  The  International  Ocean¬ 
ographic  Tables  for  the  solubility  of  oxygen  in 
seawater  J.  Cons  Cons.  hit.  Explor.  Mer  36: 
295-296. 

Senders,  j.  I...  M.  Klein,  and  J.  S.  Gallagher. 
1972.  Pres  sure- volume-temperature  relation¬ 
ships  of  gases;  virial  coefficients,  p.  4-204  to  4- 
227.  In  M.  W.  Zemanskx  [ed.|.  Heat.  Am.  Inst. 
Phvsies  Handbook.  3rd  ed..  Sect.  4.  McGraw- 
Hill. 

Truesdai.e.  G.  A..  A.  L.  Downing,  and  G.  F. 
I.owden.  1955.  The  solubility  of  oxygen  in 
pure  water  and  seawater.  J.  Appl.  Chem.  5:  53- 
62. 

Vai.entiner,  S.  1927.  01>er  die  1  .osliehkeit  der 
Edelgase  in  Wasser.  Z.  Phvs.  42  :  253-264. 

Weiss,  R.  F.  1970.  The  solubility  of  nitrogen,  oxy¬ 
gen  and  argon  in  water  and  seawater.  Deep-Sea 
Res.  17:  721-735. 

Winkler,  I..  W.  1889.  Die  Liislichkeit  des  Sauer- 
stoffs  in  Wasser.  Ber.  Dtsch.  Chem.  Ges.  22: 
1764-1774. 

- .  1891.  Die  I .osl ichkeit  der  Case  in  Wasser. 

Ber.  Dtsch.  Chem.  Ges.  24:  3602-3610. 

Submitted:  21  August  1979 
Accepted:  28  February  1980 


